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ABSTRACT 

Context. Characterising stellar and circumstellar properties of embedded young stellar objects (YSOs) is mandatory for understanding 
the early stages of the stellar evolution. This task requires the combination of both spectroscopy and photometry, covering the widest 
possible wavelength range, to disentangle the various protostellar components and activities. 

Aims. As part of the POISSON project (Protostellar Optical-Infrared Spectral Survey On NTT), we present a multi-wavelength spec- 
troscopic and photometric investigation of embedded YSOs in L 1641, aimed to derive the stellar parameters and evolutionary stages 
and to infer their accretion properties. 

Methods. Our multi- wavelength database includes low-resolution optical-IR spectra from the NTT and Spitzer (0.6-40 /jm) and photo- 
metric data covering a spectral range from 0.4 to 1 100 /jm, which allow us to construct the YSOs spectral energy distributions (SEDs) 
and to infer the main stellar parameters (visual extinction, spectral type, accretion, stellar, bolometric luminosity, mass accretion and 
ejection rates). 

Results. The NTT optical-NIR spectra are rich in emission lines, which are mostly associated with YSO accretion, ejection, and 
chromospheric activities. A few emission lines, prominent ice (H2O and CO2), and amorphous silicate absorption features have also 
been detected in the Spitzer spectra. The SED analysis allows us to group our 27 YSOs into nine Class I, eleven Flat, and seven 
Class II objects. However, on the basis of the derived stellar properties, only six Class I YSOs have an age of ~10^yr, while the 
others are older (5xlO'-10* yr), and, among the Flat sources, three out of eleven are more evolved objects (5x10^-10' yr), indicating 
that geometrical effects can significantly modify the SED shapes. Inferred mass accretion rates (Mace) show a wide range of values 
(3.6x10"' to 1.2x10"' Mq yr"'), which reflects the age spread observed in our sample well. Average values of mass accretion rates, 
extinction, and spectral indices decrease with the YSO class. The youngest YSOs have the highest Mace, whereas the oldest YSOs do 
not show any detectable jet activity in either images and spectra. Apart from the outbursting source #25 and, marginally, #20, none 
of the remaining YSOs is accretion-dominated (Lace > L,). We also observe a clear correlation among the YSO Mace, M,, and age. 
For YSOs with / > lO' yr and 0.4 Mq < M, < 1.2 Mq, a relationship between Mace and t (Mace f"' ") has been inferred, consistent 
with mass accretion evolution in viscous disc models and indicating that the mass accretion decay is slower than previously assumed. 
Finally, our results suggest that episodic outbursts are required for Class I YSOs to reach typical classical T Tauri stars stellar masses. 

Key words. Stars: formation - Stars: evolution - Infrared: stars - Accretion, accretion disks - Surveys 



1. Introduction 

Our working model of low-mass star formation arises from 
a combination of the empirical classification of the YSO 
SEDs (Lada & Wilking 1984; Lada 1987) with a theoretical pic- 
ture of YSO formation, which involves the collapse of an iso- 



* Partially based on observations collected at the European Southern 
Observatory La Silla, Chile, 082.C-0264(A), 082.C-0264(B) 



lated rotating dense core, which then forms an accreting pro- 
tostellar core and a disc (Adams & Shu 1986; Adams et al. 
1987). The empirical evolutionary sequence goes from Class 
to III objects. Class YSOs are the youngest sources (~10"^ yr), 
while Class III (the so-called weak T Tauri stars) are the old- 
est ones (10^ yr). The classification of Class 1, II, and III ob- 
jects is based on the slope of the SED between 2 and 20//m 
{ai-iQiim), defined as a = dLog(AF ,\) I dLog{A) (Lada 1987). 
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On the other hand. Class YSOs, which are usually not visible 
at these wavelengths, are defined as having LsmmlLimi >0.5%, 
where Lj,,™, is measured longward of 350/im (Andre et al. 1993, 
2000), and they have more than 50% of their mass in the sur- 
rounding envelope. The youngest protostars (0, I) are thus ex- 
tremely embedded, and are characterised by steeply rising SEDs 
from near to far-IR, mostly or entirely coming from the emis- 
sion of their surrounding envelope. According to this picture, 
most (in the case of Class I YSOs) or all (in the case of Class 
YSOs) the YSO luminosity is believed to come from accretion 
{Ltoi ~ L„cc ~ GMtMacclR*) through a circumstellar disc. Part 
of the accreted material (~10%) is ejected by means of powerful 
coUimated jets. The more evolved Class II and III sources in- 
stead have smaller IR excesses, and their SEDs can be modelled 
well by pre-main sequence photospheres surrounded by circum- 
stellar discs (i. e., classical T Tauri stars - CTTs). The luminos- 
ity of these older objects mostly originates in their stellar pho- 
tosphere, rather than from accretion processes. Thus, accretion 
and ejection activities are strongly reduced or even absent in the 
latest stages. 

Our understanding of the YSO early evolutionary stages 
(Class and I), mostly relies on SED analysis or, indirectly, on 
studies of their jets and outflows. On the other hand, we have a 
clearer picture of the latest stages, for which the physical prop- 
erties of the accreting protostars can be well studied and char- 
acterised. Recent theoretical studies on SEDs (Whitney et al. 
2003a,b; Robitaille et al. 2007) have shown that the SED anal- 
ysis alone may not be sufficient to disentangle the YSO evolu- 
tionary stage, and NIR spectroscopy is needed to characterise 
the embedded accreting protostellar core. For example, geomet- 
rical effects may produce SED misclassifications, i. e. old ob- 
jects observed edge-on may show Class I shapes, or young YSOs 
face-on may appear older, showing Class II SED shapes. Thus, 
a correct classification can be only obtained from the combined 
analysis of the YSO SED and the characterisation of the embed- 
ded stellar object and its activity. Quantitative information on 
the various phenomena characterising the environment of young 
stars can be derived from NIR spectroscopic studies, which al- 
low us to observe embedded objects, and investigate processes 
occurring in regions spatially unresolved, like the accretion fun- 
nel flows, and the ejection of jets from the adjacent inner disc, 
by using features related to the different emitting regions. 

These considerations have recently triggered a series of ob- 
servational studies, aimed at deriving the stellar physical prop- 
erties and evolutionary status of the embedded YSOs, in partic- 
ular of the Class I YSOs, which are visible at IR wavelengths. 
As a result, these studies have indicated that only a fraction of 
Class I sources is composed of highly accreting objects (see, e. g. 
White & Hillenbrand 2004; Nisini et al. 2005a; Doppmann et al. 
2005; Antoniucci et al. 2008), which show evidence of jet activ- 
ity. Several objects, classified as Class I, have indeed mass ac- 
cretion/ejection rates similar to those of Class II. Some of them 
are mis-classified, whereas others appear to be bona-fide young 
embedded objects (White et al. 2007). Indeed, the so far studied 
sample is still limited, biased, and mostly confined to a few se- 
lected star-forming regions (i. e. mostly Taurus, see, e. g.. White 
et al. 2007; Beck 2007; Prato et al. 2009). 

In this framework, we have undertaken a combined opti- 
cal/IR unbiased spectroscopic survey on a flux-limited sample 
of selected Class I/II sources, located in six different nearby 
clouds (namely Cha I & II, L1641, Serpens, Lupus, Vela, 
and Corona Australis; see, Antoniucci et al. 2011, hereafter 
paper-1), using the EFOSC2/SOFI instruments on the ESO-NTT 
(POISSON: Protostellar Objects IR-optical Spectral Survey On 



NTT). In paper-1 we presented the results of the Cha I & II 
regions, comparing La^c determinations from the different trac- 
ers, and discussing the reliability and consistency of the different 
empirical relationships considered. In this paper we report the 
survey results on L 1641, largely complemented by archive and 
literature data, which allow us to characterise the studied YSOs. 

At a distance of 450 pc, the Lynds 1641 molecular cloud 
(L 1641) is part of the Orion GMC complex (for a complete re- 
view see, Allen & Davis 2008). Southward of the ONC, L 1641 
extends from NW to SE for ~2.5°, and contains hundreds of 
young stellar objects, ranging from high- to low-mass YSOs. 
Moreover, the L 1641 cloud has been producing stars for nearly 
30Myr (Allen 1995), thus it harbours a very heterogeneous sam- 
ple of YSOs, which span from extremely active and young to old 
and quiet objects, making it the perfect candidate to study the 
different YSO evolutionary stages in an unbiased fashion. Our 
observations were performed on 27 embedded YSOs selected in 
L 1641 on the basis of the brightness and SED spectral index (a). 
Our multi-wavelength database includes low-resolution optical- 
IR spectra from the A^TT and Spitzer (0.6-40 jum), as well as pho- 
tometric data covering the spectral range from 0.4 to 1100//m, 
which allow us to construct the YSO SEDs and to characterise 
the object parameters (visual extinctions, spectral types, accre- 
tion and bolometric luminosities, mass accretion and ejection 
rates). 

This paper is organised as follows. Section 2 describes the 
selection criteria for our YSO sample. In Section 3 we define our 
observations, data reduction, and the collected literature data. In 
Section 4 we report on the results obtained from our photome- 
try and spectroscopy, we describe the detected spectral features 
and characterise the YSOs. In Section 5 we discuss the accretion 
properties of the sample, as well as the origin of the observed 
mass accretion evolution. Finally, our conclusions are drawn in 
Section 6. 

2. Sample definition 

Our sample of YSOs in L 1641 was selected from Chen & 
Tokunaga (1994) on the basis of the SED spectral index (a) 
and NIR brightness of each source, because of the instrumental 
sensitivity. In particular, we chose those objects showing a >- 
0.4 and K' <12mag, implying that most of our targets are, in 
principle. Class I and flat YSO candidates. According to te lat- 
est Spitzer surveys, this brightness constraint limits our study to 
about 15% of the entire embedded population in L 1641 (see, 
Allen & Davis 2008, and references therein). 

It is worth noting that the a classification in Chen & 
Tokunaga (1994) is based on K' and IRAS 25 fim photometry, 
and, as it will be shown in Sect. 3.3 and Sect. 4.4.2, often dif- 
fers from the 2MASS/Spitzer classification obtained in this pa- 
per (Sect. 4.4.2). This is mostly due to the low spatial resolution 
of the IRAS beam (l'x5' at 25 pm) that is not enough to prop- 
erly resolve different sources in crowded regions such as L 1641. 
As a result, the 27 YSOs are grouped into nine Class I, eleven 
Flat, and seven Class II objects. Additionally, some targets iden- 
tified as single sources in Chen & Tokunaga (1994) turned out 
to be double YSOs after analysing the Spitzer images (namely 
[CTF93]146, [CTF93]216, [CTF93]237, and [CTF93]245B). 

The targets, along with their ID, SI MB AD name, and co- 
ordinates (J2000.0), are listed in Table 1 (Columns 1-5). The 
sources have been named following the nomenclature of Chen 
et al. (1993) with the exception of those identified here as dou- 
ble YSOs, which have been labelled as " target- 1" , when coinci- 
dent or closest to the literature coordinates, and " target-2" , the 
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farthest. The large separations (>40") exclude that these sources 
are physically bound, with the exception of sources [CTF93]237, 
separated less than 5" (~2250 AU at 450 pc). Another excep- 
tion might be [CTF93]216-1 and [CTF93]216-2. They might be 
part of a wide binary system (35" ~15 700AU, at d=450pc), 
and both sources show two precessing jets in the Spitzer images 
(see also Fig. B.l in Appendix B). However such a wide system 
would not produce precessing jets, thus it is more hkely that both 
sources are binary systems. 



3. Observations and data reduction 

Our multi-wavelength database on L 1641 targets is made of: i) 
low-resolution optical and NIR spectroscopic observations from 
the ESO/NTT, to infer most of the YSO main stellar parame- 
ters (spectral type, accretion, mass accretion and ejection rates); 
ii) archival Spitzer-IRS spectra, mostly used to derive the visual 
extinction; Hi) Spitzer-IRAC and Spitzer-MIPS images from the 
Spitzer Heritage Archive^, and literature photometric data, cov- 
ering a spectral range from 0.4 to llOO/zm, which allow us to 
construct the SEDs. 

3. 1 . Optical and NIR spectroscopy 

Our optical and near-infrared spectroscopic observations 
were obtained at the ESO New Technology Telescope with 
EFOSC2 (Buzzoni et al. 1984) and Sofl (Moorwood et al. 1998), 
respectively. The observations were carried out over a short pe- 
riod of time (10-13 and 14-16 February 2009, for the NIR and 
optical spectroscopy, respectively), thus possible YSO variabil- 
ity should not significantly affect the optical and NIR segments 
of our spectra. The observational settings (spectral resolution 
and slit width) were chosen to obtain, as far as possible, homo- 
geneous spectra, aiming to apply a combined optical/NIR anal- 
ysis. The total integration time on source (/,) in both optical and 
NIR spectra differs from target to target, depending on the source 
brightness, which was chosen to obtain a signal to noise ratio 
>200 on the continuum and detect emission lines, which are 1/20 
of the continuum, with at least S/N =10. No preferential position 
angles for the slits were chosen. 

The EFOSC2 spectra were taken with the grism N.16 and a 
(J.'l slit width, covering a wavelength range from ~0.6 to 1 fim, 
and with a spectral resolution of R~700. /, ranges from 30 s (tar- 
gets with y=l 1 mag) up to 1800 s (targets with y=17 mag). Only 
those targets having a V band magnitude <17 were observed 
(namely 14 out of 27 objects), and they are indicated in our ob- 
servation log of Table 1 (column 6) with the label 'O'. 

The Sofl NIR spectra were taken with an AB offset scheme. 
We used a 0'.'6 slit (R~900) for both the blue grism (ranging 
from 0.95 to 1.64yum, and labelled as 'B' in column 6 of Tab. 1) 
and the red grism (ranging from 1.51-2.5//m, and labelled as 
'R' in column6 of Tab. 1). Only those targets with J<15mag 
were observed with the blue grism (namely 16 out of 27 tar- 
gets), whereas all the targets were observed with the red grism. 
/, ranges from 20 s (objects with J=8mag) up to 1200 s (ob- 
jects J>13mag) in the blue grism, and from 20 s (targets with 
K=5 mag) up to 1800 s (K>12 mag) in the red grism. 

Additionally, telluric and spectro-photometric standards 
were observed to correct for the atmospheric spectral response 
and to flux-calibrate the spectra, respectively. In particular 
HD 289002 (B3 spectral type) was chosen for the optical spec- 



' http://sha.ipac.caltech.edu 



tra, whereas an F8V type star (Hip 23930) was selected for the 
NIR. 

The data reduction was done using standard IRAF^ tasks. 
Each spectrum segment was flat fielded, sky subtracted and cor- 
rected for the curvature derived from long-slit spectroscopy, 
while atmospheric features were removed by dividing the spec- 
tra by the telluric standard star The wavelength calibrations were 
performed using a helium-argon lamp in the optical and a xenon 
lamp in the infrared. The resulting spectra were originally flux 
calibrated by means of the observed standard stars. We thus ob- 
tained for each target up to three individually calibrated spectra 
(optical, blue and red grism). However, due to variable seeing 
conditions during our observations (from (^.'6 to r.'3) and Sofl 
instrumental problems (the source did not remain correctly cen- 
tred on the slit during the nodding cycle), often the stellar cali- 
brated continua do not properly match in the three spectral seg- 
ments. Thus, we performed aperture photometry on the Ks band 
Sofl acquisition images (see Sect. 3.3), to get an absolute flux- 
calibration for the red grism spectra, and match to these the stel- 
lar continua of the two remaining segments, finally obtaining a 
single flux-calibrated spectrum from 0.6 to 2.5 fim. Uncertainties 
on the Ks photometry range from ~5 to 10% of the total flux. 

3.2. Archival Spitzer spectroscopy 

From the Spitzer Heritage Archive we retrieved 20 spectra (the 
spectrum of [CTF93]50 - #2 - has already been presented by 
Fischer et al. 2010), which were labelled as 'S' in column 6 of 
Tab. 1. [CTF93]237-1 and -2 are not spectrally resolved, thus 
the spectra are blended. No additional data are available in the 
archive for the remaining objects in the Spitzer archive. 

All the data were taken with the Spitzer Infrared 
Spectrograph (IRS; Houck et al. 2004) under a variety of pro- 
grams (namely 30706, 30859, and 50374), and observed be- 
tween November 2006 and November 2008. All the IRS spectra 
were obtained with a combination of short-low (SL2/SL1, 5.2- 
8.7/7.4-14.5;um, R = 64-128), and long-low (LL, 14.0-39.0//m, 
R - 64-128) modules (i.e. ranging from 5 to 39 /^m). Source 
#9 was observed with the LL module alone. The total expo- 
sure time on source ranges from four to eight minutes. All 
the spectra were observed in IRS staring mode and extracted 
from the Spitzer Science Center SI 8.7.0 pipeline basic calibrated 
data (BCD) using the Spitzer IRS Custom Extraction (SPICE) 
software. This process includes bad-pixel correction, optimal 
point-spread function (PSF) aperture extraction, defringing, or- 
der matching, wavelength and flux calibration. The extracted 
spectra from the eight nod positions were averaged together 
(four for the SL2/SL1 and LL modules, respectively), obtaining 
a final flux-calibrated spectrum, which ranges from ~5 to 39 i-im. 

3.3. Imaging: Sofl, Archival Spitzer imaging 

Acquisition images in the Ks band were taken with Sofl before 
each NIR spectroscopic observation, with an /, of one minute 
each. The data were reduced using IRAF packages and applying 
standard procedures for sky subtraction, dome flat-fielding and 
bad pixel removal. 

Aperture photometry was performed on the resulting images 
using the task phot in IRAF. The images were, then, flux cal- 



- IRAF (Image Reduction and Analysis Facility) is distributed by 
the National Optical Astronomy Observatories, which are operated 
by AURA, Inc., cooperative agreement with the National Science 
Foundation. 
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ID 


Target name 


Other Name 


('• 


(J2000.0) 

HI 


5(J2000.0) 

r ' ") 


S pec trosc opy" 


1 


[CHS2001] 13811 


2MASS-J05360665-0632171 


05 


36 


06.65 


-06 


32 


17.1 


0,B,R,S 


2 


[CTF93]50 


NAME HH 147MMS 


05 


36 


25.13 


-06 


44 


41.8 


0,B,R,S 


3 


[CTF93]47 


V* V380 Ori 


05 


36 


25.43 


-06 


42 


57.7 


0,B,R 


4 


[CTF93]32 


V* V846 Ori 


05 


36 


41.34 


-06 


34 


00.3 


0,B,R,S 


5 


[CTF93]72 


IRAS 05350-0700 


05 


37 


24.47 


-06 


58 


32.9 


0,B,R,S 


6 


[CTF93]83 


HH 43 IRSl 


05 


38 


07.44 


-07 


08 


30.0 


R,S 


7 


[CTF93]62 


V* V1787 Ori 


05 


38 


09.23 


-06 


49 


15.9 


0,B,R,S 


8 


[CTF93]79 


V* V883 Ori 


05 


38 


18.10 


-07 


02 


25.9 


0,B,R 


9 


[CTF93]99 


Haro 4-254 


05 


38 


52.36 


-07 


21 


09.5 


0,B,R,S 


10 


[CTF93]87 


2MASS-J05390536-071 1052 


05 


39 


05.36 


-07 


11 


05.2 


0,B,R,S 


11 


[CTF93]104 


HBC 176 


05 


39 


22.32 


-07 


26 


44.5 


0,B,R 


12 


Meag31 


2M AS S- J0540 1 379-0732 1 59 


05 


40 


13.79 


-07 


32 


15.9 


B,R 


13 


[CTF93] 146-2 


2M AS S- J0540 1 494-0748485 


05 


40 


14.94 


-07 


48 


48.5 


R,S 


14 


[CTF93] 146-1 


IRAS 05378-0750 


05 


40 


17.80 


-07 


48 


25.8 


R,S 


15 


[CTF93]211 


IRAS 05379-0815 


05 


40 


19.40 


-08 


14 


16.3 


0,B,R,S 


16 


[CTF93]187 


V* V1791 Ori 


05 


40 


37.36 


-08 


04 


03.0 


0,B,R,S 


17 


[CTF93J168 


IRAS 05389-0759 


05 


41 


22.13 


-07 


58 


03.0 


R,S 


18 


[CTF93]191 


V* DL Ori 


05 


41 


25.34 


-08 


05 


54.7 


0,B,R,S 


19 


[CTF93]246A 


2MASS-J05413005-0840092 


05 


41 


30.05 


-08 


40 


09.2 


R,S 


20 


[CTF93]186 


IRAS 05391-0805 


05 


41 


30.22 


-08 


03 


41.4 


R,S 


21 


[CTF93]246B 


2MASS-J05413033-0840177 


05 


41 


30.33 


-08 


40 


17.7 


R,S 


22 


[CTF93]237-2 


2MASS-J05413418-0835273 


05 


41 


34.18 


-08 


35 


27.3 


R,S 


23 


[CTF93]237-1 


IRAS 05391-0836 


05 


41 


34.78 


-08 


35 


23.0 


R,S 


24 


[CTF93]216-1 


IRAS 05403-0818 


05 


42 


47.06 


-08 


17 


06.9 


R,S 


25 


[CTF93]216-2 


2M AS S- J05424706-08 1 7069 


05 


42 


48.48 


-08 


16 


34.7 


R,S 


26 


[CTF93]245B-2 


[BDB2003]G212.98-19.15 


05 


42 


50.48 


-08 


38 


29.2 


0,B,R 


27 


[CTF93]245B-1 


IRAS 05404-0841 


05 


42 


50.50 


-08 


39 


57.6 


B,R 



Notes. Labels O, B, R, and S indicate EFOSC2 optical spectra, SofI blue-, and red-grisms, and Spitzer spectra, respectively. See text for details. 



ibrated using several field stars as 'standards' along with their 
Two-Micron All Sky Survey (2MASS, Skrutskie et al. 2006) 
band photometry. An aperture correction was estimated on a few 
isolated stars using the task mkapfile in IRAF with a formal er- 
ror of 0.02-0.03 mag. As a result, we get uncertainties for the 
YSO photometry in the range of ~0. 05-0. 2 mag, depending on 
the number of standards and the YSO brightness. 

In addition, Spitzer basic calibrated data (BCD) have been 
obtained from the Spitzer Heritage Archive, reduced with the 
S 14.0 pipeline. 

They consist of several individual frames, which map 
the entire L 1641 cloud, observed with the Infrared Array 
Camera (IRAC; Fazio et al. 2004) in four channels (at 3.6, 4.5, 
5.8, and S.O/^m) (covering a region of ~1.04deg^) and with the 
Multiband Imaging Photometer for Spitzer (MIPS; Rieke et al. 
2004) (at 24, 70, and 160 yum) (-1.5 deg^). The IRAC images 
were performed in High-Dynamic Range mode with integration 
times of 0.4 and 10.4 s (Spitzer program ID 43), and have been 
presented in previous papers (Megeath et al. 2005; Fang et al. 
2009). The MIPS images (Spitzer program ID 47) have effec- 
tive integration times of 80, 40, and 8 s at 24, 70, and 160yum, 
respectively, and have been presented by Fang et al. (2009). 

A mosaicking of the individual frames in the seven different 
bands was performed using the Spitzer MOPEX (MOsaicker 
and Point source Extractor) tool (Makovoz et al. 2005), which 
performs background matching of individual data frames, mo- 
saics the individual frames, and rejects outliers. IRAC and MIPS 
24 jum photometry for our targets was reported by Megeath et al. 
(2005) and Fang et al. (2009), and kindly provided by these au- 
thors. 



We thus performed photometry only on the MIPS 70 and 
160 jum reduced mosaics. We manually examined the latter for 
target detection and used MOPEX/ APEX for the source extrac- 
tion. We detect 26 out of 27 targets at 70 yum, and 1 1 out of 27 in 
the 160/im map. 

As indicated in the MIPS data handbook, we used apertures 
of radius 16" and 32" for the 70jum and 160/im photometry 
respectively, a sky annulus with inner and outer radii of 18" 
and 39" at 70yum, and 64" and 128" at 160/im and an aper- 
ture correction factor of 2.07 and 1 .98 for the 70 pm and 160 jum 
photometry, respectively. Uncertainties were obtained from the 
pixel-to-pixel noise in the sky annulus. This error is added in 
quadrature to the photometric uncertainty. Source identification 
and extraction were particularly difficult for the 160//m map, 
due to the extreme brightness and saturation of the 160//m dif- 
fuse emission in several parts of the map, which translated into 
very large photometric errors or into a non-detection of the tar- 
gets. 

3.4. Additional photometry from the literature 

Additional NIR photometry in the J, H, and Ks bands was re- 
trieved from 2MASS (observed in November 1998), with 10 cr 
detection limits of 16.2, 15.3, and 14.6 mag, respectively. 

Additional photometric data for the sample were re- 
trieved from 'the Naval Observatory Merged Astrometric 
Dataset' (NOMAD; Zacharias et al. 2004) in the optical which 
includes USNO optical photometry (B, V, R, I bands). More op- 
tical data were taken from the Sloan Digital Sky Survey (SDSS 
York et al. 2000), which includes measurements at u'g'r'i'z' 
bands at 3540, 4760, 6290, 7690, and 9250A. 
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Table 2. 2MASS and SofI band magnitudes of the observed targets 



ID 


K^ilMASS) 


K (So fl) 


AAr,(5 - 2M)° 










1 


9.68±0.02 


9.3+0.1 


-0.38 


2 


8.21±0.02 


8.2+0.1 


-0.01 


3 


5.95+0.02 


6.10+0.05 


0.15 


4 


9.52+0.03 


9.3+0.1 


-0.22 


5 


9.85+0.02 


9.8+0.1 


-0.05 


6 


11.79±0.07 


11.6+0.2 


-0.19 


7 


7.98+0.02 


7.6+0.06 


-0.38 


8 


5.15+0.01 


5.7+0.05 


0.55 


9 


8.03+0.03 


8.00+0.06 


-0.03 


10 


10.54+0.02 


10.5+0.1 


-0.04 


11 


8.17+0.03 


8.1+0.05 


-0.06 


12 


11.1+0.02 


11.2+0.2 


0.1 


13 


13.47+0.06 


13.5+0.2 


0.03 


14 


10.49+0.02 


10.3+0.1 


-0.19 


15 


10.22+0.02 


10.2+0.1 


-0.02 


16 


7.93+0.05 


7.90+0.05 


-0.03 


17 


10.53+0.02 


10.6+0.1 


0.07 


18 


9.29+0.02 


9.20+0.07 


-0.09 


19 


12.27+0.03 


12.4+0.2 


0.13 


20 


10.56+0.02 


10.1+0.1 


-0.41 


21 


11.86+0.03 


11.8+0.2 


-0.05 


22 


11.52+0.02 


11.6+0.2 


0.08 


23 


10.57+0.04 


10.1+0.1 


-0.47 


24 


11.76+0.03 


11.8+0.1 


0.04 


25 


11.06+0.03 


7.98+0.05 


-3.08 


26 


9.09 +0.02 


9.2+0.1 


0.11 


27 


10.36+0.02 


10.5+0.1 


0.14 
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Fig. 1. Spitzer/MIPS 70yum map of L 1641. Crosses indicate the sample 
targets. 



K.iSofI)- KA2MASS). 



Finally sub-MM photometry was retrieved from Di 
Francesco et al. (2008) (SCUBA/JCMT, at 450 and 850;um), 
and from Dent et al. (1998) (UKT14/JCMT, at 350, 800, 850, 
and 1100 //m). 

4. Results 

4.1. Imaging 

Table 2 reports the 2MASS, and SofI Ks band photometry for 
the sample, as well as the YSO variability between the SofI and 
2MASS observations (corresponding to about ten years). From 
a quick inspection of Column 4 in Tab. 2, we first note the large 
Ks variation (~3.1mag) of source #25 (namely [CTF93]216- 
2), whose outburst was recently reported by our team (Caratti o 
Garatti et al. 201 1). Five more objects (namely # 1, 7, 8, 20, 22) 
show some degree of variability (0.3< AKs <0.5mag), whereas 
the remaining twenty objects do not indicate any, inside their 
photometric error bars. 

Tables A.l and A. 2 (reported in Appendix A) present the 
complete collected photometry for our targets from 0.4 to 
11 00 yum. We do not consider and report the available IRAS mea- 
surements in this table, due to the crowding of the fields at mid- 
IR wavelengths, and the availability of the Spitzer data. 

In Fig. 1 we report the Spitzer-MIPS map at 70jum, indicat- 
ing the location of the studied objects. 

4.2. Optical and NIR spectra 

Optical and NIR continuum-normalised spectra of the sample 
are shown in Figures 2 and 3. The most prominent features ob- 



served in the spectra are emission lines, which have been indi- 
cated in the figures. The fluxes of the main lines, used for the 
analysis in this paper, are reported in Tables 3 and 4, along with 
the statistics on the detections. A complete list of all the de- 
tected lines in each spectrum, including fluxes, full width half 
maximum (FWHM), and equivalent widths (Wa), is given in 
Appendix A (Tables A.3-A.29). The equivalent widths and line 
fluxes were calculated by integrating across the line, after sub- 
tracting the continuum, which was estimated by interpolating be- 
tween two line-free adjacent regions. 

The detected lines are circumstellar features and are mostly 
associated with YSO accretion activity or inner winds, as, e. g., 
H I, Ca II, He i (see, e. g., Muzerolle et al. 1998; Natta et al. 2004; 
Edwards et al. 2003; Antoniucci et al. 2008), or ejection activity, 
as, e.g., [Oi], [Sii], [Fell], H2 (see, e.g., Hartigan et al. 1995; 
Nisini et al. 2002; Caratti o Garatti et al. 2006). A few YSOs, 
namely #2, #3, # 1 1, and #18, also show permitted ionic emis- 
sion lines (e.g., Oi, Fez, Feii, Nai, Mgii, Ci), which are char- 
acteristic of active young stars (see, e. g., Hamann & Persson 
1992a,b; Kelly et al. 1994; Hernandez et al. 2004), and usuaUy 
associated with chromospheric activity. In particular, the Fei, 
and Fell emission lines have been observed in several young 
eruptive stars during their outburst activity (e.g. V1647 Ori, 
see, e.g., Fedele et al. 2007), or in YSOs showing high mass- 
accretion rates (e. g., Rossi et al. 1999). 

The most commonly detected features are from H i recom- 
bination lines, Ho' in the optical spectra (11 out of 14, i. e. 79% 
detection rate), Paschen lines between 0.8 and 1.3 fim (e. g. Pa/3 
has a 69% detection rate - 10 out of 16), and Brackett lines be- 
tween 1.5 and 2.2 jum (Bry has a 96% of detection rate - 26 out 
of 27). Other prominent bright emission lines are the Ca 11 triplet 
between 0.85 and 0.87 fim (detected in 57% of the optical spec- 
tra), and the He: line at 1.08/^m (56% detection rate), which 
notably displays a clear P Cygni profile in four out of nine detec- 
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tions (namely source #9, 10, 15, 18), usually indicative of inner 
winds in accretion disc systems (see, e. g., Edwards et al. 2003), 
with radial velocities larger than 300 km s"' . Source # 2, namely 
[CTF93]50, only shows the blue-shifted Hei line in absorption, 
and no indication of emission. Connelley & Greene (2010) sug- 
gest that this particular feature may be indicative of FU Orionis 
like stars, and might have a disc wind origin. According to Fang 
et al. (2009), the optical spectrum of [CTF93]50 (their source 
# 105) resembles those of FU Ori-like YSOs, however the typ- 
ical absorption features of the CO band-head lines longward of 
2.29jum, also characteristic of FU-Ori objects, are not detected 
in our NIR spectrum. 

Lines tracing the jet activity have also been identified in our 
sample, but with a detection rate lower than that of the accretion 
lines. Namely the [O i] line at 0.63 fim has been identified in 43% 
of the optical spectra, the [S ii] doublet at 0.68 yum in 36%, the 
[Fe II] 1 .64 ;um and H2 2. 12//m lines in 19% and 32% of the NIR 
spectra, respectively. In total, jet-tracer emission lines have been 
observed in 14 out of 27 spectra, i. e. in 52% of the sample. 

Our low resolution spectra mostly lack absorption lines on 
the photospheric continuum, also indicating the presence of veil- 
ing. Source #7 and # 16 are an exception, and clearly show Hi 
absorption features at optical and NIR wavelengths, indicating 
both early spectral types and lower veiling'', likely indicating 
a more evolved stage. On the other hand, broad-band molecu- 
lar absorption features, typical of cool objects, are clearly visi- 
ble in both optical (VO and TiO bands) and NIR (H2O bands) 
spectra. CO overtone bands have been detected in our spectra, 
both in absorption (sources #8, #11, #18, #24) and emission 
(# 14 and #25). The CO bandheads in emission are usually be- 
lieved to originate in the inner gaseous disc (see e. g. Carr et al. 
1993; Antoniucci et al. 2008), and have been observed in YSOs 
that display strong jets and Herbig-Haro (HH) objects (see, e. g., 
Davis et al. 201 1), whereas the CO bandheads in absorption may 
originate in cooler regions, like the YSO envelope (Davies et al. 
2010) or the photosphere of cool YSOs (see e. g. Rayner et al. 
2009). 



For those sources where absorption features such as Hi (spectral 
types F or earlier), Nai and/or Cai (spectral types K or M) have been 
detected, we usually get < 1. For the remaining sources, we can infer 
lower limits for the veiling, on the basis of non-detection of Nai or Cai 
lines in the K band (see, e. g., Antoniucci et al. 2008). Assuming MO as 
average spectral type and typical intrinsic equivalent widths ~ 3. 3 A, we 
get r<. > 1. 
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Fig. 2. Continuum-normalised spectra, EFOSC and SofI Z band. Labels indicate the most prominent spectral features. 
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Table 3. Main detected emission lines: accretion tracers. 



Source 


Hff 




Can 0.854/0.866 /jin 


He I 1.083/jm 




Pad 




Pay 




Pa/3 




Bry 


1 


ID 


(F± AF)10-" 


EW 


(F± AF)10"" 


EW 


(F± AF)10-"' 


EW 


(F± AF)10-''' 


EW 


(F ± A F)10-''' 


EW 


(F± AF)10-" 


EW 


(F±AF)10-'*' 


EW 




erg s"' cm"^ 


A 


erg s"^ cm"^ 


A 


erg s"^ cm"" 


A 


erg s"^ cm"" 


A 


erg s"^ cm"" 


A 


erg s"' cm"" 


A 


erg cm"- 


A 


1 


13±1 


-11.0 










0.8±0.3'' 


-2.0 


1.1±0.3 


-2.8 


1.8±0.4 


-4.0 


3.0±0.5 


-5.0 


2 


0.99±0.04'' 


-7 


1.8±0.2/1.6±0.2 


-6.4/-5.7 






1.1±0.2 


-2.5 


2.0±0.4 


-3.3 


4.5±0.4 


-4.4 


7.4±0.8 


-3.3 


3 


1330±1 


-73.2 


344±5/304±5 


-26.6/-23.7 


95±8 


-5.1 


152±8 


-8 


283±8 


-16.2 


386±8 


-23.2 


181±8 


-13.0 


4 


82.3±0.2 


143.1 


2.3±0.4/3.1±0.4 


-4/-5.5 


9±2 


-8.5 


6±2 


-5.7 


8 ±2 


-8.1 


10±1 


-9.2 


3.3+0.8 


-4.5 


5 


0.13±0.05" 




0.41±0.03/0.38±0.03 


-19.0/-17.5 


0.33±0.08 


-3.7 






0.89±O.O8 


-9.5 


2.8±0.1 


-12.7 


3.8±0.4 


-7.6 


6 


























0.4±0.1 


-6.5 


7 
8 


18±1 


-6.3 






















12±3 
140±10 


-3.1 
-3.6 


9 


13.4±0.2 


-17 






3.8±0.8'' 


-1.9 










42±1 


-1.6 


10±1 


-3.8 


10 


3.05±0.03 


-96.1 


1.2±0. 1/1.1 ±0.1 


-18.6/-14.1 


0.89±0.08'' 


-4.5 


0.5±0.1 


-3.3 


1.3±0.2 


-6.1 


2.2±0.1 


-8.5 


2.2±0.4 


-8.9 


11 


69.9±0.1 


-88.2 


26.4±0.4/23.0±0.4 


-27.9/-24.4 


28.8±0.7 


-17.5 


19±2 


-14.7 


27.0±0.8 


-16.1 


45±1 


-21.2 


30±1 


-13.1 


12 


























0.6±0.1 


-4.3 


13 


























0.2±0.1° 




14 


























3.0±0.3 


-11.1 


15 


5.18±0.05 


-95.6 


1.35±0.4/I.06±0.6 


-21.1/-16.5 


l.l±0.l'' 


-8.7 


0.5±0.1 


-5.5 


1.2±0.1 


-9.5 


2.0±0.2 


-10.0 


2.5±0.4 


-7.2 


16 


























7±1 


-2.5 


17 


























1.8±0.3 


-6.9 


18 


403.0±0.5 


102.9 


42.2±0.6/35.8±0.6 


-44/-36 


21.5±1'' 


-9.7 


29±1 


-13.1 


38±1 


-17.6 


45±1 


-25.2 


14.1±0.5 


-16.5 


19 


























0.3+0.1 


-9.8 


20 


























2.7±0.4 


-7.1 


21 


























0.7±0.1 


-7.3 


22 


























<0.1 




23 


























3.8±0.9 


-5.0 


24 


























1.74±0.08 


-12.4 


25 


























10±2 


-3.6 


26 










0.40±0.09 


-3.9 


0.2±0.1" 


-3.8 


0.4±0.1 


-4.0 


1.5±0.2 


-5.5 


4.1 ±0.6 


-5.2 


27 


























1.6±0.8° 


-11.1 


Detections 


1 1/14 - 79% 


8/14-57% 


9/16 - 56% 


9/16-56% 


10/16-62% 


11/16-69% 


26/27 - 96% 



Table 4. Main detected emission lines: ejection tracers. 



Som"ce 


[0 1] 0.630//m 




[S ii] 0.672/3 iim 




[Feu] 1.644//m 


H2 2.122//m 


1 


ID 


(F± AFjlO"" 


EW 


(F± AF)10-''' 


EW 


(F± AFjlO"" 


EW 


(F± AF)10-"' 


EW 




erg s"' cm"" 


A 


erg s-2.3"' cm"" 


A 


erg s"' cm"" 


A 


erg s"' cm"" 


A 


2 


0.77±0.02 


-7.1 


0.33±0.03/0.55±0.03 


-2.1/-3.4 


2.1 ±0.8° 


-1.1 






4 


1.3±0.1 


-2.7 














5 










2.9±0.2 


-1.9 


1.4±0.4 


-2.8 


6 














0.2±0.1° 


-5.4 


9 














3±1 


-1.3 


10 


0.26±0.02'' 


17.2 


0.07±0.0370.09±0.03 


-2.6/-3.1 






0.9±0.4° 


-3.5 


11 


2.0±0.1 


-3.1 


0.4±0.1 


-0.5 










12 














0.7±0.1 


-4.3 


14 










0.57±0.07 


-6.6 






15 


0.54±0.05 


12.6 


0.26±0.04/0.32±0.04 


-4.5/-5.5 






0.8±0.2 


-2.5 


18 


9.1 ±0.4 


-2.6 


1.5±0.4/3.0±0.4 


-0.4/-0.8 


4±1 


-3.2 






22 














0.2±0.1'' 


-2.8 


24 










2.48±0.05 


-45.8 


0.7 ±0.06 


-4.9 


Detections 


6/14-43% 


5/14 - 36% 


5/27-19% 


8/27 - 30% 



Notes. S/N ratio below 3(t. Line with P Cygni profile. 
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4.3. Spitzer IRS spectra 

Our IRS Spitzer spectra are shown in Figure 4. The flux den- 
sities for the Y-axes have been scaled and shifted for a better 
display. Moreover, the three panels have a different Y-axis range 
to properly fit the different steepness of the spectra. Such a vari- 
ety already shows that the SED spectral index of our YSOs has a 
spread larger than originally expected (see Sect. 2), i. e. > -1, 
instead of o- > -0.4, as we will see later in Sect. 4.4.2. 

There are a few emission features detected in the Spitzer 
spectra, mostly H2 pure rotational lines (at 17.0 and 28.2 //m) 
(detected in 50%, i.e. 10/20), [Feii] (at 17.9, 22.9, 24.5, and 
35.8 jum) (detected in 40%, i. e. 8/20), [Si 11] and [Si iii] (at 34.8 
and 35.3 jum, and at 25.6 fim) (detected in 50% and 55%, respec- 
tively, i.e. 10/20 and 11/20), and [Siii] (at 33.5yum) (detected 
in 65%, i. e. 13/20). There are also a few lines which we tenta- 
tively assign to [Fe i] (at 24, 28.5, and 36.5 fim) (detected in 30%, 
i. e. 6/20). These lines are quite common in Mid-IR YSO spec- 
tra (see, e. g., Forbrich et al. 2010; Lahuis et al. 2010; Baldovin- 
Saavedra et al. 2011), and mostly originate from both dissocia- 
tive (e.g. [Sm], [Sin], [Sim], [Feii], [Fei]) and non-dissociative 
(H2) shocks along the YSOs jets (Lahuis et al. 2010), whereas a 
few lines may also originate from the discs, as H2, [Si 11], [Si iii], 
[Fei] (see, e.g., Lahuis et al. 2007). A complete list of all the 
detected lines is given in Appendix A (Tables A.3-A.29). 

Several YSOs in our sample show prominent ice absorp- 
tion features at 5-8 //m (H2O) and IS.lfim (CO2) (79% of the 
sample) as well as the amorphous silicate absorption feature at 
9.7 fim (observed in all the IRS Spitzer spectra). This last feature 
is detected in absorption in all sources except #7 and #18, where 
it is observed in emission, usually indicating more evolved ob- 
jects (see, e. g., Watson et al. 2004). These broadband features 
are usually connected to dust and ice mantels on dust in the in- 
terstellar medium and in the YSO envelope, and their optical 
depths have been used in several works to compute extinction 
towards the protostellar photospheres (see, e.g., Alexander et al. 
2003;Chiar et al. 2007). 

4.4. Derived Stellar parameters 
4.4.1. Reddening 

To obtain correct estimates of important stellar parameters, such 
as the accretion luminosity {Lace), bolometric and stellar lumi- 
nosities {Lhoi and L»), mass accretion and mass ejection rates 
{Mace, and Mom), the observed line fluxes and photometry must 
be properly dereddened. Thus, deriving an accurate value of the 
visual extinction (Av) towards the stellar photosphere and its cir- 
cumstellar region is a fundamental task. There are several inde- 
pendent methods to compute Ay, and they are successfully ap- 
plied to CTTs, that are usually less embedded. Determining the 
extinction towards more embedded protostars is particularly dif- 
ficult, however, due to the high extinction and/or scattered light 
from the outflow cavities and reflection nebulae surrounding the 
YSOs (see, e. g., Whitney et al. 2003b; Beck 2007; ConneUey & 
Greene 2010). 

In this work, we have combined our multi-wavelength spec- 
troscopy and photometry to obtain independent estimates of the 
visual extinction (see Table 5), namely: 

;) Line ratios of transitions arising from the same upper 
level can be used to evaluate Ay. Indeed the observed ratio de- 
pends only on the differential extinction, if the emission is opti- 
cally thin, and once the theoretical value, which depends on the 
Einstein coefficients and frequencies of the transitions, is known. 



L1641 




[5.8]-[8.0] 

Fig. 5. Top panel: JH vs. HK diagram for our L1641 targets. Dwarf MS 
and classical T Tauri loci (Meyer et al. 1997) are indicated by a contin- 
uous black line and a dashed magenta line, respectively. Black dotted 
dashed lines show the reddening vectors. Bottom panel: IRAC colour- 
colour diagram for the L1641 targets. The black dotted dashed rectangle 
delineates the domain of Class 11 YSOs (Allen et al. 2004), two paral- 
lel lines indicate the locus of reddened Class 11 sources (Megeath et al. 
2004), and an average reddening vector for Av=30mag is plotted, as 
well. 



We use the Bry to Pa^ ratio, assuming that the observed 
emission arises from optically thin gas. Results are reported in 
column 2 of Tab. 5. [Feii] line ratios are also used, in partic- 
ular 1. 644/1. 257//m (source #2, #5, and #18), 1.644/1. 32 l//m 
(source #5) or 1.71/1.60yL(m (source #24), adopting the transi- 
tion probabilities of Nussbaumer & Storey (1988). Results are 
reported in column 3 of Tab. 5. 

ii) The optical depth of the 9.7/im amorphous silicate ab- 
sorption feature (rg 7) and optical depth of the CO2 ice feature at 
15.2yum (T15.2) can also be used to estimate the extinction. 

To measure the optical depths, the following steps were fol- 
lowed. First we determined the stellar continuum by fitting a 
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Fig. 4. Spitzer IRS spectra. Flux densities have been shifted and normalised by a constant value for a better display. Labels indicate the detected 
spectral features. 



low-order polynomial on a log-log scale to the spectral segments 
between 5 and 35 jum, which are not affected by absorption fea- 
tures (see e. g. Boogert et al. 2008). Then, we applied to derive 
T,i, the following relationship (see, e.g., Alexander et al. 2003; 
Chiar et al. 2007): 

= -ln(F,oi,rceW/FcontW) (1) 

where F source and F continuum are the flux densities of the source 
and of the continuum fit at A, respectively. T^(C02),ce can be con- 
verted to a column density, using the following relationship (see, 
e. g., Alexander et al. 2003); 

N(C02) = (2) 

peak 

where J T^dA is the equivalent width, Apgak is the wavelength 
of the feature peak optical depth, and A is the ice band strength 
(1.1 X 10'^ cmmolecule"') assuming pure CO2 ice (Gerakines 
et al. 1995; Alexander et al. 2003). 

The Ay value for the silicates is obtained from Rieke & 
Lebofsky (1985) as Ay- Tg 7X 18, and from the CO2 column den- 
sity using the relationships of Bergin et al. (2005). Results from 
the CO2 and silicate measurements are reported in Columns 4 
and 5 of Tab. 5, respectively. 

Hi) Colour-colour diagrams from J-, H-, and /T-band 
2MASS photometry. Fig. 5 (top panel) shows the J - H \s 
//-/T colours of our targets, along with the dwarf-MS and CTTs 
loci (Meyer et al. 1997), and the reddening vectors. Many of the 
sources are located along the reddening band extending from the 
CTTs locus, while there are a few objects positioned on the right 
side of the strip, indicating less evolved objects. Arrows on these 
objects indicate lower limits, because their /-band magnitude is 
an upper limit or has a S/N ratio <3. Derived Ay values are re- 
ported in column 6 of Tab. 5. 

iv) Self-consistent method (SCM, see Sect. 4.4.5 and 
Antoniucci et al. 2008, for a detailed description). This method 
relies on the fact that there are two ways of determining the 
source stellar luminosity L,, both depending on the actual 
extinction value: i) from the bolometric luminosity Lhoi (see 
Sect. 4.4.2) and the accretion luminosity Lace derived from the 
de-reddened Pa/? (and/or Bry) flux (Sect. 4.4.3), assuming that 
L, - Lboi - Lace', ii) from the observed K band magnitude of 
the source, considering the distance, the spectral type, and an es- 
timate of the veiling in the relevant band (see Sect. 4.4.4). An 



average veiling value of r^: = 1 for the whole sample is assumed. 
Because veiling dependence on magnitude is logarithmic, while 
extinction dependence is linear, the computation is less sensitive 
to veiling variations. For example, assuming = 5, instead of 
one, would increase Ak of ~0.48, i. e. ~5 mag in the V band. 

The best estimate of the extinction will be the one for which 
we obtain the same value of L, from both computations. 

For comparison, additional extinction estimates retrieved 
from the literature (Strom et al. 1989; Fang et al. 2009; 
Connelley & Greene 2010) are also reported in column 8 of 
Tab. 5. 

The inferred Ay values range from to ~30mag, although 
the different methods adopted here usually produce a wide range 
of results for each object. Indeed this is a well known problem 
for Class I sources (see e. g. Beck 2007; Prato et al. 2009; Davis 
et al. 201 1), and it is due to several reasons, as e. g.: ;) the consid- 
ered lines trace different regions of the YSO (e. g., the circum- 
stellar region or the jet, the H i or the [Fe 11] lines, respectively), 
possibly with different extinctions; //) scattering is not taken into 
account by any of the above methods, although those based on 
the [Fe 11] line ratios or on MIR optical depths should not be sig- 
nificantly affected; ///) a standard ISM law (Rieke & Lebofsky 
1985) has been assumed to correct for the differential extinction, 
without considering that it can vary, depending on the size and 
properties of the grains in the disc, envelope, or jet (see, e. g. 
Cardelli et al. 1989). 

Despite these limitations. Tab. 5 shows a relatively good 
agreement among the different methods. In particular, the cross 
correlation coefficient between different sets of Ay ranges be- 
tween r = 0.95 (Si vs SCM) with an rms of ~3 mag, to r = 0.66 
(CC vs SCM) with an rms of ~6mag. The average value is 
r - 0.86. We then adopt an average Ay value for each object 
(column 9 of Tab. 5), after discarding those that deviate from the 
average by more than one sigma. 

4.4.2. SEDs and Lboi 

Both photometric and spectroscopic data of our sources are in- 
corporated in the SEDs presented in Figure 6, where both target 
IDs and names are also reported. In general, there is an excellent 
agreement between observed optical-IR spectra and photometric 
points, independently obtained from images. However, there is a 
marginal discrepancy between the two datasets for a few objects, 
especially at short wavelengths (see, e. g., source #6, # 13, # 19), 
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Table 5. Visual extinction values of the sample. 
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which is likely due to scattering effects, and/or intrinsic YSO 
variability. Mismatches between Spitzer spectra and photometry 
are also observed in a few sources (see, e.g., #5, #10, #14), 
and might be due to variability, as already observed in the liter- 
ature (see, e. g., Morales-Calderon et al. 2009). Special mention 
is needed of source #25 (or [CTF93]216-2), which underwent a 
strong outburst between 2007 and 2008 (Caratti o Garatti et al. 
201 1), increasing in brightness by ~4.6, 4.0, 3.8, and 1.9 mag in 
the /, H, Ks bands and at 24 jum. The panel in Figure 6 shows 
both pre- and outburst SEDs. Photometric pre-outburst and out- 
burst values are reported in Tab. A. 1 and in Caratti o Garatti et al. 
(201 1), respectively. 

The YSO Class is determined from the SED slope (a, or 
spectral index), which is calculated over all the available data 
points between 2.2 and 24yum. The spectral index value is ob- 
tained from a linear fit to the logarithms, taking into account 
the uncertainties in the flux measurements. Depending on a, 
our sources were classified as Class I (a > 0.3), Flat spectrum 
(0.3 < a < 0.3), or Class II (1.6 < a < 0.3). a values of the 
sample range from 1.60 to -1.04, including 9 Class 1, 11 Flat, 
and 7 Class II YSOs. Notably, as already observed for other out- 
bursting sources (see, e.g., Kospal et al. 2011, and references 
therein), source # 25 changed its colours and its SED slope (from 
0.69 to 0.25, pre- and outburst phase, respectively). Despite its 
actual a values, in this paper we consider it as a Class I source. 

The source ID, their a values and corresponding classifica- 
tions are reported in Columns 1, 2, and 3 of Table 7, respec- 
tively. As akeady mentioned in Sect. 2, the a values (and thus the 
YSO classification) listed in Chen & Tokunaga (1994) often dif- 
fer from ours, in some cases substantially as, e. g., for source #9 
(-0.45 vs 1.1, our and their value, respectively) or for source #1 1 



(-0.49 vs 1), indicating that IRAS fluxes should not be used for 
those objects located in crowded regions, as the resulting classi- 
fication is not reliable. 

The total bolometric luminosity {Lboi) of each source was 
derived by integrating the observed SED (see Figure 6), which 
includes data ranging from optical to mm wavelengths. The cal- 
culation was performed starting from the first available SED data 
point (which varies from source to source, see also Table A.l), 
using straight line interpolation (in the Log{X) — LogiAF^) plot) 
between all the SED points, also including additional points 
from the observed spectra. A final correction at the longest wave- 
length was applied assuming that F^i decreases as A^^. A distance 
value of 450pc was assumed for the entire sample (see, e.g., 
Allen & Davis 2008). Li,oi inferred values for the sample range 
from ~0.3 to 188 Lq, and are listed in column 4 of Table 7. For 
the outburst source #25, both pre- and outburst Lhoi and a values 
are reported. 

It worths noting that no correction for geometrical effects 
has been applied to the infeiTed bolometric luminosities, because 
system inclinations are not well defined. On the other hand, SED 
models from Whitney et al. (2003b) indicate that, due to geomet- 
rical effects, the observed bolometric luminosities of embedded 
objects may differ from their intrinsic luminosities up to a factor 
of two between pole-on and edge-on systems. Another source of 
uncertainty for Lhoi may be given by the foreground extinction, 
which cannot be coiTectly estimated, because the infeiTed Ay 
values of Sect. 4.4. 1 do not discriminate between foreground and 
circumstellar contribution. However, extinction maps can pro- 
vide us with upper limits for the foreground extinction. From 
Rowles & Froebrich (2009) maps, we infer an average Ay value 
of ~6 mag towards our targets, which would modify Lhoi values 
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Fig. 6. Target SEDs and spectra. Photometric data with error bars are in black, spectroscopic data (ESO-NTT and Spitzer) are in red. Target IDs 
and names are also reported. - to be continued 
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Fig. 6. Target SEDs and spectra. Photometric data with error bars are in black, spectroscopic data (ESO-NTT and Spitzer) are in red. Target IDs 
and names are also reported. - continued 



between 5 and 10%, because foreground extinction only affects 
the shortest wavelengths of the SED. As a consequence, esti- 
mates of Lhoi and other stellar parameters (i. e., L,, M,, age, 
and Mace), that will be later derived from it, might also be af- 
fected by such uncertainties. 



4.4.3. Accretion luminosity 

The accretion luminosity is the main indicator of YSO accretion 
activity. Many different techniques exist in the literature to derive 
it (see, e. g., Calvet et al. 2000), from the measurement of the U- 
band excess luminosity (GuUbring et al. 1998) to the luminosity 
of optical and infrared emission lines (as, e.g.. Ha, Can, Pa/3, 
Bry; e. g., Muzerolle et al. 1998; Calvet et al. 2004; Natta et al. 
2004, 2006), which are thought to be mainly produced in the 
magnetospheric accretion flow. The inferred relationships have 
been successfully tested for a wide range of YSO masses, from 
substellar-mass to intermediate-mass YSOs (Calvet et al. 2004; 
Natta et al. 2004; Garcia Lopez et al. 2006), but have been rarely 
tested simultaneously on the same YSO sample (Muzerolle et al. 
1998; Rigliaco et al. 2011), because of the wide wavelength 
spread between the considered lines and the intrinsic YSO vari- 
ability, which require simultaneous coverage of the optical and 
infrared spectral wavelength range. 

In paper-1 we analysed various empirical \mt-Lacc relation- 
ships from five different tracers, namely [O i] at 6300A, Wa, Ca ii 
at 8542A, Pa/3, and Bry, critically discussing the various deter- 
minations in the light of the source properties. As a result, we 
showed that the Bry and Pa/? lines give the smallest dispersion 
of Lace over the entire range of L«, whereas the other tracers, 
especially the Ha and [O i] lines, provide much more scattered 
Lace results, not expected for the homogeneous sample of targets 
observed. 

Here, we use our 0.6-2.5yL(m spectra to derive Lace simulta- 
neously from the Can, Pa/?, Bry line luminosities. 



As in paper-1. Lace from these lines is derived from the fol- 
lowing empirical relationships (Herczeg & Hillenbrand 2008; 
Calvet et al. 2000, 2004): 



Log{Lace/LQ) = 1-02 X Log(Lcam542lLQ) + 2.5 
LogiLacclLo) = 1.03 X LogiLpaplLo) + 2.80 
Log(LacelLQ) = 0.90 x Log{LBry/LQ) + 2.90 



(3) 
(4) 
(5) 



Line luminosities were obtained after de-reddening the ob- 
served fluxes by the adopted Ay (Table 5), using the standard 
Rieke & Lebofsky (1985) extinction law, and assuming a dis- 
tance of 450 pc. Results are reported in Table 6, where source ID, 
Lflfc(Bry), LacdPaf)), and Lflcc(Caiio.854) are listed. The derived 
values range from ~0.1 to ~60Lo. We assigned an upper limit 
to source #23, computed from the upper limit for Bry emission. 

In Figure 7 we compare the obtained results, plot- 
ting Log{Laee(Paf5)) vs Log(Lacc{Bry)) (black circles), and 
LogiLacciCau)) vs Log{Lacc{Bry)) (red triangles). The dashed 
line indicates the equivalent locus. In the first case. Pa/? vs Bry, 
there is a very good agreement with a correlation coefficient 
r equal to 0.92, whereas the correlation of the second dataset 
(Can vs Bry) gives a worst match with r - 0.8, as can also be 
noted from the larger scatter in the plot. Such a difference can- 
not be attributed to a wrong extinction estimate, with the pos- 
sible exception of source #11 (see also Tab. 6), otherwise, due 
to differential extinction, we would observe the data points be- 
low the dotted line with the circles positioned above the triangles 
(overestimated extinction), or the opposite (underestimated ex- 
tinction). It is worth noting that our low resolution spectroscopy 
does not allow us to de-blend the Can (0.854/im) and the Pa 15 
(0.855/im) lines, thus the observed Can flux at 0.854yL(m might 
be overestimated up to ~20%. As a consequence, we prefer to 
discard the values from the Ca n lines. For each source we adopt 
an average Lace value, obtained by averaging the accretion lumi- 
nosities inferred from the Bry and Pa^ Unes. Results are reported 
in column 7 of Table 7. 
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Table 6. Accretion luminosity estimates from Bry, Pa/?, and 
Can 0.854 //m lines. 
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Fig. 7. Comparison between the Lace values derived from different 
lines: Pa/3 vs Bry (black circles), and Can vs Bry (red open triangles). 
Correlation coefficients are also reported. 



4.4.4. Stellar spectral types 

We obtain the spectral types (SpTs) of our targets by combin- 
ing various methods. The derived SpTs along with the method 
used and the effective temperatures {T^s) are reported in Table 7 
(Columns 6 and 7, respectively). The adopted T^s is converted 



from the spectral type adopting the relations given by Kenyon 
& Hartmann (1995) (for SpTs earlier than MO), and by Luhman 
et al. (2003) (for MO or later SpTs). 

For the 14 sources with optical spectra (see Sect. 3.1 and 
Tab. 1) we employ the spectral classification code SPTclass de- 
veloped by Hernandez et al. (2004)"*. This code uses empiri- 
cal relationships between the equivalent widths (EWs) of many 
atomic/molecular absorption/emission lines and Teif . As a result, 
this automated code can classify optical spectra to a precision 
of about one subtype. Additionally, spectral classifications for 
many of these optical sources (9 out of 14) already exist in the lit- 
erature (Strom et al. 1989; Fang et al. 2009; Connelley & Greene 
2010). There is a good match (~1 subtype) with our results, ex- 
cept for the possible FU-Ori object [CTF93]50 (our source #2, 
see also Sect. 4.2), which shows a later spectral type (K6+1) 
with respect to previous measurements from Strom et al. (1989) 
(G5-G7, source KMS 31) and Fang et al. (2009) (K3-4.5, source 
#105 in their work). Indeed these differences are not surprising, 
because FU-Ori sources may show large spectral type variabili- 
ties (see e. g. Hartmann & Kenyon 1996). 

On the other hand, the thirteen remaining NIR spectra (i. e. 
sources #6, 12, 13, 14, 17, 19, 20, 21, 22, 23, 24, 25, 27) show 
a steeply rising flat continuum that is almost featureless or with 
absorption water bands (in the J, H and K bands, between 1.3- 
1.5 /^m, and 1.7-2.1/^m). These bands arise at T^ff <3800K (MO 
or later spectral types), and thus are typical of M type objects. 
Thus, these thirteen spectral types have been identified following 
the prescription given by Gorlova et al. (2010), i. e. dividing the 
spectra in two groups, which either show or do not the absorption 
water bands in the NIR. 

H2O broad band features are clearly detected in eight out of 
thirteen spectra (namely #12, 13, 19, 21, 22, 23, 24, 25). The 
depth of these bands increases with decreasing temperatures, 
thus the values of their EWs, once corrected for the proper Ay, 
can be related to the spectral subtype (see, e.g., Kleinmann & 
Hall 1986; Greene & Lada 1996; Luhman et al. 2003; Gorlova 
et al. 2010). Usually the 1.3-1.5/(m absorption feature is de- 
tected for Al/ <20mag, and the EW of the 1.7-2.1 //m feature de- 
creases by about a factor of 2, as the visual extinction increases 
from to 30 mag {Av-25 mag is the maximum extinction value 
in our sample, see also Tab. 7, column 5). Along with the absorp- 
tion water bands, source #12 (Meag31) also shows VO bands in 
absorption at ~1.15/(m (typical of M6 or later SpTs; see, e.g.. 
Gushing et al. 2005). On the other hand the lack of strong narrow 
band features, mostly due to the veiling and the low resolution of 
our spectra, prevents us from an accurate spectral classification. 
For these objects, the uncertainties on the SpT range from one to 
two subtypes, i. e. up to ~200-300K. 

The five remaining spectra (namely sources #6, 14, 17, 20, 
27) are featureless with steeply rising SEDs, and without any 
sign of overturn in the H or K band, which means either SpT 
later than MO, or SpT earlier or equal to MO and Av>30-40 mag 
(see, e.g., Greene & Lada 1996; Gorlova et al. 2010). Our Ay 
estimates exclude the latter hypothesis, thus we safely classify 
them as K spectral types. One of them ([CTF93]245B-1, i.e. 
source #27) was already classified as K6 by Strom et al. (1989), 
thus we keep this classification. The lack of features in the re- 
maining four sources precludes any detailed spectral analysis. 
Therefore, we use both our spectral and photometric data (see 
also Fig. 6) to model the four YSOs with the SED fitting tool 
from Robitaille et al. (2007) and provide some tighter constraints 
on the underlying stellar objects. Indeed, the SED fitting tool 



http://www.astro.lsa.umich.edu/ hernandj/SPTclass/sptclass.html 
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Fig. 8. Temperatures inferred from SED fitting. For each YSO the stel- 
lar temperature distribution (hashed histogram), normalised to its max- 
imum, is plotted over the entire grid of models (in grey). 



alone cannot confidently predict the stellar temperature in em- 
bedded sources. Even knowing the L» value, there is still degen- 
eracy between and values. This degeneracy can be par- 
tially removed, given L«, Lboi, extinction, and distance. In fact a 
substantial increase (decrease) of (500 K or more) produces 
a decrease (increase) of 7?,, generating an older (younger) YSO 
and substantially modifying the SED shape. 

Thus, as further constraints for the fitting tool, distance, ex- 
tinction, and bolometric luminosity of each source were used (as 
listed in Tab. 7). An additional constraint is given by the YSO 
spectral type (K), and thus by the temperature range. Then, for 
each object, we obtained a grid of possible models, listed as a 
function of their values. We then select the best 50 models 
for each YSO, plot the output stellar temperature distribution, 
and infer the most likely temperature. The distribution spread 
gives us a rough estimate of the error. Our results are reported in 
Fig. 8, where, for each object, the stellar-temperature distribu- 
tion (hashed histogram), normalised to its maximum, is plotted 
over the entire grid of models (in grey). 

The derived spectral types for the entire sample are reported 
in Table 7 (column 6), and range from B7 to M7.5. Most of the 
objects (~74%) are low-mass YSOs with SpT between K5 and 
M7.5, four (-15%) are intermediate mass YSOs (F7 to K4), and 
three (-11%) are Herbig Ae/Be stars (A5 to B7). 

4.4.5. H-R diagram and stellar parameters 

To characterise our YSOs it is necessary to derive the stellar 
parameters (L,, M,, Rt). L, of each source is inferred from 
L, = Lhoi - Lace, i- e. assuming that the observed Lboi is sim- 
ply the sum of accretion and stellar luminosities. The resulting 
L, values are reported in Table 7 (column 8). Assuming that the 
YSO SED is well modelled by our data, i. e. that the derived Lboi 
has a relatively small uncertainty, the error on L, is due to the 
uncertainty of the Lace estimates. 

To test the consistency of these results we derive L* from 
the stellar bolometric magnitude (see, e. g., Gorlova et al. 2010), 
obtained from the dereddened Ij and J magnitudes (Table A.l, 
Columns 5 and 6, respectively) for 22 out of 27 and 26 out of 27 
objects of the sample: 

Log(LJLQ) = 1.86- 0.4(7 -0.482Av + BC/- 5/0^(450) H-5)(6) 
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Fig. 9. L, values derived from the Ij and / band photometry (red tri- 
angles and black dots, respectively) versus L, values obtained from 
L, = Li,„i - L„cc- The dashed line marks the equivalue locus. 
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Fig. 10. HR diagram for the whole sample. The evolutionary models 
of Siess et al. (2000) are shown with isochrones at lO**, 10^, 10*, and 
5x10* yr (red dashed lines), and mass tracks from 0.1 to 3Mq (black 
dashed lines). 



LogiLJLQ) = 1. 86 -0.4(/-0.265Av + BC7- 5/0^(450) -h5)(7) 

where the bolometric correction in Ij and J [BCj and BCj, 
respectively), for each SpT of Table 7, are inferred from the 
dwarf colours in Kenyon & Hartmann (1995). 

Figure 9 compares in logarithmic scale L, values obtained 
from I J and J magnitudes (red triangles and black dots, respec- 
tively) against L, = Lboi -Lace values of our Table 7. The dashed 
line marks the equivalue locus. Both samples show some scatter- 
ing around the equivalue locus, but, as expected, values derived 
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Fig. 11. Spectral index (a) versus age. The two dashed lines delineate 
the Class 1-Flat (a = 0.3) and Flat-Class II (a = -0.3) boundaries. The 
object IDs are also reported. 

from the J band have a larger scatter and tend to be located above 
the dashed line. This is because eq. 7 does not take into account 
the infrared excess in the J band, which is relevant in young stel- 
lar objects (see e. g. Cieza et al. 2005), thus L, is overestimated. 
On the other hand, the Ij band should be much less affected by 
the infrared excess, and thus the derived L, value should be more 
reliable, if the extinction has been correctly estimated. 

Our Lt and T^ff estimates can be plotted on a HR diagram 
to infer M,, R^, and age of each sampled source. We adopt the 
evolutionary tracks from Siess et al. (2000), with a metallicity 
of Z = 0.02, and we use the on-line tool^ to infer these results 
and produce the HR diagram. The HR diagram for the sample 
is shown in Figure 10, and the inferred M«, and age values 
are reported in Table 7 (Columns 9, 10 , and 11, respectively). 
For self-consistency, isochrones and ages down to -lO'^yr are 
reported in both Figure 10 and Column 1 1 of Table 7, however, 
we stress that age estimates below ~10'' yr are not reliable. Thus, 
the reader should consider these values more like a qualitative in- 
dication of stellar youth, rather than a real estimate. Moreover, it 
is worth noting that adopting different sets of evolutionary tracks 
would provide different values for the derived stellar masses and 
ages, up to a factor of 2-3 (see, e.g., Hillenbrand et al. 2008; 
Spezzi et al. 2008; Fang et al. 2009), where the largest discrep- 
ancies are in the age estimates of objects with age <1 Myr For 
example, for low-mass objects, the morphology of Siess et al. 
(2000) tracks is quite similar to those of Baraffe et al. (1998), 
thus there is good agreement between ages, masses and radii 
given by the two models (~20-40 %). On the other hand, for ex- 
ample, tracks from D' Antona & Mazzitelli (1997) give, on aver- 
age for our sample, lower mass and younger age estimates than 
Siess et al. (2000). As a consequence, these quantities and those 
later inferred, like, e.g., the mass accretion rate in Sect. 4.4.6, 
can differ up to a factor of 2-3, depending on the model. 

According to the Siess et al. (2000) model, the sampled stel- 
lar masses range from 0.4 to 3.4 M0, with an average and me- 

' http://www-astro.ulb.ac.be/ siess/database.html 



dian values of -1.1 and 0.9 Mq, respectively. The average and 
median age of the sample is about 2 and 1 Myr, respectively, 
with values ranging from <10^ to 10^ yr In Figure 1 1 the mea- 
sured a values, derived in Sect. 4.4.2, are plotted against stellar 
ages, indicating that, on average, a decreases with stellar age. 
However, there is no straightforward correlation between spec- 
tral index and age. Six out of nine Class I YSOs have an age 
of ~10^, while the remaining three (namely #6, 19, and 20) ap- 
pear to be older (5x10^ to 10^ yr). The estimated age for the 
flat spectrum sources ranges from ~10^ to ~10^ yr, with an aver- 
age value of 2 Myr, whereas the age of Class II sources spans 
from -10* to -lO^yr, with an average value of 4 Myr. This 
indicates that the standard SED classification (a computed be- 
tween 2.2 and 24 //m) might not properly reflect YSO age, i. e. 
that the SED slope might not be a good indicator of the stellar 
age. Indeed, geometrical effects may play a role in modelling 
the YSO SED (see, e. g., Whitney et al. 2003b; Robitaille et al. 
2007), altering the slope and generating misclassifications. This 
is particularly clear when the SEDs show a double peak (at opti- 
cal/NIR and MIR wavelengths), often indicative of edge-on discs 
or transitional discs (see, e.g., Cieza et al. 2007; Merm et al. 
2010; Williams & Cieza 201 1, see also Fig. 6). 

Finally, we also note that two massive sources, namely #7 
and #16, are ~10^yr old. These are the only objects showing 
strong H i absorption lines in the optical/NIR sampled spectra 
(see Sect. 4.2 and Fig. 2 and 3). Since they seem to be older 
than L 1641 (-5x10'' yr; eg., Allen & Davis 2008), this might 
indicate that they are not part of this star forming region. 

4.4.6. Mass accretion rates 

Mass accretion rates (Mace) for the whole sample can be derived 
once the accretion luminosities and the stellar parameters have 
been inferred. Since the accretion luminosity corresponds to the 
energy released by the accreting matter onto the YSO, assuming 
that the free fall starts at the co-rotational radius (7?,„), i. e. at 
~5R, (Gullbring et al. 1998), L^cc is then: 

L,,, - G/W./W«,,( 1 - IRin)IR, (8) 

and thus Mgcc is given by: 

Macc^L„,,^\.25R,IGM, (9) 

The derived Mgcc values are reported in Table 7 (Column 14). 
These values span four orders of magnitude, ranging from 
3.6x10"' to 1.2x10"^ Mgyr"', with the highest Mgcc given by 
the low-mass outbursting source #25 ([CTF93]216-2; Caratti o 
Garatti et al. 2011), and the lowest value (source #23) being 
an upper limit. Error estimates are particularly difficult, because 
they depend on both observational and theoretical uncertainties, 
namely on Lgcc, R*, and M,. Therefore an error bar up to one 
order of magnitude can be expected. 

Mass accretion estimates for five targets of the sample have 
been previously reported in literature (see also notes in Table 7). 
Fang et al. (2009) characterised four sources, namely #1, #2, 
#4, and #18, by means of optical spectroscopy. Estimates for 
two sources (# 1 and #18) are in good agreement with our re- 
sults, whereas the others (sources #2, and #4), differ of about 
one order of magnitude, possibly because of the different Ay and 
Lboi estimates. On the other hand, source #2 has also been anal- 
ysed by Fischer et al. (2010), who modelled photometric and 
spectroscopic data (Spitzer and Herschet), infer an Mgcc value 
very similar to ours. Finally, Hillenbrand et al. (1992) derive 
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source # 3 (e. g. V* V380 Ori) stellar properties from optical and 
IRAS photometry, inferring a mass accretion rate of ~ one order 
of magnitude higher than ours, hkely because they obtain larger 
stellar mass and luminosity. 

4.4.7. Jet detection and mass ejection rates 

Accretion and ejection activity are complementary and intercon- 
nected processes of the stellar birth. In particular, part of the ac- 
creting material is ejected by the YSO in the form of a bipolar 
jet, which, in turn, sweeps up the circumstellar and interstellar 
medium, producing a molecular outflow (see, e. g., Reipurth & 
Bally 2001). Strong accretion is usually accompanied by strong 
ejection processes, thus the detection and study of jets/outflows 
give us insights into YSO activity. The YSO outflow-activity is 
evinced from typical jet-tracing lines in the spectra (see also 
Sect. 4.2 and Tab. 4), narrow-band imaging centred on these 
lines, and/or CO maps tracing wider outflows. 

With these ideas in mind we inspected our Spitzer images in 
search of possible signatures of jets from the sampled sources. 
Indeed, the IRAC bands contain both molecular hydrogen and 
ionic lines, which may be shock-excited in protostellar outflows. 
In particular, band 2 (4.5 yum) contains bright molecular hydro- 
gen lines and can be used to detect shock spots (knots) (see, 
e. g., Peterson et al. 201 1). To better identify these features in the 
IRAC/Spitzer mosaics, we thus constructed IRAC three-colour 
images, using 3.6//m, 4.5 //m, and 8/zm (i. e. bands 1, 2, and 4, 
in blue green, and red, respectively), and identified the knots by 
means of their colours and morphologies. We detect jets from 
eleven sources (namely sources #2, 3, 5, 6, 8, 10, 11, 15, 18,24, 
25), nine of them were already observed by Davis et al. (2009) in 
their H2 survey of L 1641, whereas sources #24, and #25 were 
not covered by that survey. A faint H2 emission was discovered 
by Connelley et al. (2007) around #24. The two jets from #24, 
and # 25 detected in the IRAC/Spitzer mosaics are described and 
shown in Appendix B. Spitzer images usually reveal bright jets 
emitting in the mid-IR, but, can fail to detect faint jets, only vis- 
ible in narrow-band images. Therefore, we also searched the lit- 
erature for indications of further jets/outflows from our YSOs. 
In Column 16 of Table 7 we indicate the presence of jets and 
outflows driven by our sources, as reported by Allen & Davis 
(2008); Davis et al. (2009); Connelley et al. (2007) (jets), by 
Morgan et al. (1991); Dent et al. (1998) (CO outflows), and by 
this paper, whereas in Column 17 we indicate the presence of 
jet-tracing lines in our optical/NIR spectra. 

Additionally, to compare accretion and ejection activity in a 
statistical way, we also infer a crude estimate of the mass ejec- 
tion rates from the different lines tracing jets (namely [O i], [S 11], 
[Fe 11], and H2 lines reported in Tab. 4), and observed in our opti- 
cal/NIR spectra (see Sect. 4.2). As these lines are optically thin, 
their luminosity gives us an estimate of the total mass (M). Thus 
the mass ejection rate {Mom) can be inferred if the tangential ve- 
locity (^t) and jet extension {It) are known, i.e. Mom - Mvt/lt- 
Indeed our low resolution spectra do not allow us to measure the 
radial velocity, neither do we know the inclination angle of the 
jets, thus we assume an average Vt (150kms"' and 50kms"', for 
the atomic and H2 lines, respectively; see, e.g., Hartigan et al. 
1995; White & Hillenbrand 2004; Caratti o Garatti et al. 2009). 
Moreover, the jet size is assumed equal to the measured seeing, 
projected to the L 1641 distance. This assumption rests on the 
fact that the aperture-extraction width in our spectra is defined 
by the seeing limited width of the stellar continuum, thus the jet 
is not spatially resolved, and on the fact that the spectroscopic 
absolute flux calibrations were done using the photometry. 



For the atomic species we use the following relationship, 
Moui - pmn{NuV)vtlk (e. g. Nisini et al. 2005b; Caratti o Garatti 
et al. 2009), with A^hV = Lx(/!vAi/if |||)-', where Lx is the lu- 
minosity of the element X, for the selected transition, Ai and /i 
are the radiative rate and the fractional population of the upper 
level of the transition, y is the ionisation fraction of the consid- 
ered species with a total abundance of |p| with respect to hydro- 
gen. We assume that the element X is completely ionised, and 
element abundances for Orion from Esteban et al. (1998) and 
Esteban et al. (2004). For the [Oi] (6300A) and [Sn] (6731 A) 
lines, we follow prescriptions given by Hartigan etal. (1995), us- 
ing their equations A8 and A 10, respectively. To derive the [Fe 11] 
line (1.64yL(m) intensities, for all the sources we assume an elec- 
tron density of n^-lQr' cmT^, i. e. close to the [Fe 11] critical den- 
sity. This particular value is typical of the jet base (Takami et al. 
2006; Garcia Lopez et al. 2008, 2010), and it has been obtained 
from source # 24, using the different [Fe 11] line ratios observed 
(namely 1.64/1.53, 1.64/1.60, 1.64/1.66, and 1.64/1. 68pm), and 
adopting the technique used by Nisini et al. (2002) and Takami 
et al. (2006). Fluxes of optical and NIR lines were dereddened 
using the adopted Ay values reported in Tab. 7 and the deredden- 
ing law form Rieke & Lebofsky (1985). 

In a similar way, the value of Mout(H2) can be written as 
M01AH2) - l^im.YiN'n.AvJk (e. g. Nisini et al. 2005b; Caratti o 
Garatti et al. 2009), where is the average atomic weight, 
the proton mass, A^h, the H2 column density, A the area of the H2 
knot (i.e. encompassed by the slit), the tangential velocity, and 
/t the projected length of the knot (in this case the slit width). The 
A^Ht value has been obtained from the dereddened intensity of the 
1-0 S(l) line (2. 12yum), assuming a typical shock temperature of 
2 000K. 

We infer an average Mom value for those sources with more 
than one mass-loss rate estimate. The derived Mo,„ values are 
reported in Column 15 of Table 7. Because of the several as- 
sumptions made, these are very crude estimates with at least one 
dex uncertainties, and they should not be considered accurate for 
single sources. The Mom estimated values range from 4x10"'° to 
6x10"'' M0 yr"', i. e. from one to two orders of magnitude lower 
than those of well known powerfuU jets and HH objects (e. g., 
Podio et al. 2006; Antoniucci et al. 2008; Caratti o Garatti et al. 
2009). The mass ejection rates estimated here are thus more con- 
sistent with less powerful outflows, driven by less powerful ac- 
cretors, as also observed by Hartigan et al. (1995) in Taurus or by 
White & Hillenbrand (2004) in the Taurus and Auriga regions. 

4.4.8. Statistics 

Once we have derived the main physical parameters of the YSOs 
(Tab. 7), we can then analyse our sample in a statistical way, and 
investigate whether, and how, these quantities vary depending on 
different types of subsamples. 

To this aim, we have selected six subsamples: a) three sub- 
samples were selected according to the class of the sources 
(Class I, Flat, and Class II); b) a subsample called 'Outflow 
YSOs', that includes the sampled YSOs with an outflow/jet 
signature in the images (Col. 16 of Tab. 7) and/or in the opti- 
cal/NIR spectra (Col. 17 of Tab. 7); c) a subsample called 'Jetless 
YSOs', that includes the sampled YSOs with no outflow/jet sig- 
nature; d) a subsample called 'Very young YSOs', that includes 
the sampled YSOs with age <5xlO^ yr (Col. 1 1 of Tab. 7). The 
full sample as well as the subsamples do not include the out- 
burst source #25, because of its peculiar characteristics, which 
would affect the statistics. As main observables we take into ac- 
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count Mace, Mom, MacclMout ratlo, stellar age, LacdUoi, Ay, a, 
and Mt, for which an average value is derived. The number of 
each subsample elements («) is in round brackets, and it dif- 
fers for M„„i, where values are reported according to Column 
15 of Tab. 7 Although some of these parameters have large un- 
certainties (e. g., Mout) when considering single objects, they are 
still significant for statistical purposes. Results are reported in 
Table 8. Mace, as well as Ay, a, and M» decrease with YSO class, 
whereas the average age increases, as expected. The mass ejec- 
tion rate decreases from Class 1 to Flat, while it increases in the 
Class II sample. However, this value is computed from only two 
Class II sources (which show jet lines), thus it does not represent 
a proper statistical sample. Sources with outflows have mass ac- 
cretion rates ~1 dex higher than 'jetless' YSOs, which are older, 
indicating that accretion/ejection activity decreases with time. 
The MacclM„„i ratio seems to be constant (~0.01), although, due 
to the uncertainties, this value is not very indicative. On the other 
hand, the LacdLboi ratio is slightly higher in younger and active 
sources. 
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Table 7. Stellar Parameters 



ID 




Class 




Av 


SpT 




L, 


M, 


R, 


Age 






Mace 




Jet/Outflow'' 


Jet lines'" 








(Le) 


(mag) 




(K) 


(Le) 


(Mo) 


(Rq) 


(yr) 


(L©) 




(Moyr-') 


(Moyr-') 






1 


1 OA 


TT 
11 


1 ^ 


L 


JV J ± 1 




n Q 

\j.y 


1 1 7 
1.1/ 


1 .7 J 


4 71 P-i-A 


n 4 


U.J 1 


9 1 F S 




n 


n 


Z 


u.uo 


r 


1 A ^ 


O.J 


JVD ± 1 


A7nn-i- 1 


1 ^ 


1 1 ^ 
i . i J 


A 

D.'+J 


9 1 9P-i,^ 
Z. iZJ_,+J 


1 .J 


n 1 

U. lU 


9 7F 7 
Z. /C,- / 


7T^ no 


y 


y 


acc 
J 


-\J.Hy 


TT 
11 


U 1 .T- 


9 


/\1 ±1 


^Z^UltZJU 


zL9 A 


9 AS 
Z.OO 


9 AS 
Z.Oo 


4 9AP-I.A 


1 S S 
1 0.0 


U.J 1 


A OF 7 
o.uc- / 




y 


n 


Add 




r 


^ 1 


J 


iVlu ± i 


Jo jU± 1 jU 


9 

Z.J 


U.J J 


^ 49 
J.4Z 


J. lVJi+ J 


A 
U.D 


1 Q 
\J. Ly 


1 9F 7 




n 


y 


c 

D 


n 1 n 


r 


J. / 


1 n 


i^^*" -I- 1 
j\.j ± 1 




9 

Z.O 


1 1 

1 . 1 


9 Q1 


Q 7 1 

y. 1 1 J_,+ J 


1 1 

1 . 1 


^ 
U.J 


Q ^F S 


or? no 


y 


y 


U 


U.J / 


T 
1 






i\ / ^ It J 




J. 1 


AO 


1 SA 
J. JU 


A ^AF-I-S 

u. Jum-T J 


s 

U.J 


14 

U. IH- 


S 9F S 


AV 10 


y 


y 


7 


U.J / 


TT 
11 


1 ^ 1 


u 


/\ J IE i 


5700+900 


1 1 Q 
11." 


1 .00 


1 SA 

1 .OH- 


1 OOF 4-7 


1 9 

1 .z, 


OQ 
u.u^ 


^ 7F S 

J. /li-O 




n 


n 


Q 

o 


u.uo 


p 


1 oo. J 




J57r.f + J 


1 ^000+ 1 000 


127 


3.43 


9 7Q 


9 SlF-i-A 


61.3 


0.33 


1 6F-6 

1 .Uli u 




y 




9 


-0.45 


JJ 


6 


2 


K4' •'' + 1 


4'^qO+l so 


4.4 


1.51 


3.23 


9.62E+5 


1.6 


0.27 


l.lE-7 


2E-09 


y 


y 


10 


0.28 


F 


3.0 


8.2 


MO' ' ± 1 


3850+130 


2.4 


0.55 


3.69 


4.13E+5 


0.6 


0.20 


1.3E-7 


8E-10 


y 


V 

J 


11 


-0.49 


JJ 


7.1 


2 


YSf +2 




4.4 


1.89 


2.39 


5 OlF+6 


2.7 


0.38 


l.lE-7 


2E-09 


y 


y 


12 


0.33 


J 


1.1 




]yj7 5/ + J 5 


7800+700 


0.9 


0.15 


3.71 


1.76E+4 


0.2 


0.18 


1.6E-7 


8E-10 




y 




1.60 


J 


Q 


zj 


]yjj/ +2 


^705+^00 

J / UJX JUU 


J .u 


47 


4 40 


5.01E+4 


^ 

U.J 


0.08 


Q OF-8 

7.UJ_- 








14 


0.50 


J 


7.0 


24 


K7* ±3 


3890+400 


3.9 


0.60 


4.10 


3.24E+5 


3.1 


0.44 


6.8E-7 


5E-09 


y 


y 


15 


-0.04 


F 


2.3 


8.5 


K7'' ± 1 


4060+ 1 00 


1.7 


0.73 


2.87 


7 66F+5 


0.6 


0.26 


7.5E-8 


7E-10 


y 


V 

y 




-0.26 


F 


10.2 


4 


P7<:.rf.c _,_2 


6780+160 


9.1 


1.70 


2.54 


y . J T u 


1.1 


0.1 1 


5.2E-8 








17 


0.16 


F 


4.7 


23 


K5* ±2 


4T',0+300 


2.9 


1.12 


2.94 


s Q^F+S 


1.8 


0.38 


1.5E-7 








18-^^ 


0.27 


F 


5.8 


2 


J^jr.rf.f + J 


5080+170 


3.7 


1.74 


2.14 


5.34E+6 


2.1 


0.36 


8.3E-8 


6E-09 


y 


y 


19 


0.39 


J 


1.4 


20 


MO' ±2 


-3C';n+3oo 

J ^ V/ J ' \J\J 


1.1 


0.56 


2.25 


1 oOF+6 


0.3 


0.21 


3.8E-8 






JJ 


20 


0.32 


I 


4.7 


23 


K59 +2 


4T',0+300 


2.2 


1.12 


2.31 


1 ^OF+6 


2.5 


0.53 


2.1E-7 






JJ 


91 


U.Zj 


r 


z.u 


77 5 


\/rfV -t-7 

IVlu ItZ 


^SSO-i-^00 


1 9'^ 

1 .Zj 


'^A 
U. JO 


9 

z. jy 


0. JZIi^r J 


7S 
U. / J 


"^7 
U. J / 


1 OF 7 
1 .uc- / 




n 


n 


22 


-0.14 


F 


0.32 


12 


MO' ±2 


3850+300 


0.26 


0.59 


1.09 


6.00E+6 


<0.06 


0.19 


<3.6E-9 




n 


n 


23 


0.90 


I 


7.7 


21 


MO'' ±2 


3850+300 


4.9 


0.60 


5.26 


6.01E+4 


2.8 


0.36 


7.9E-7 


4E-09 


n 


y 


24 


0.93 


I 


6.8 


25 


M2^' +2 


3560+300 


4.7 


0.4 


5.73 


1.O6E+4 


2.1 


0.31 


9.6E-7 


6E-09 


y 


y 


25 


0.69/0.25^ 


I/F 


2/22'' 


18 


M5^±1.5 


3200+200 


1.9 


0.24 


4.40 


2.13E+4 


20 


0.91 


1.2E-5 




y 


n 


26 


-0.56 


II 


2.6 


10 


K5''-'^ ±2 


4350+300 


1.8 


1.11 


2.25 


1.63E+6 


0.8 


0.31 


6.3E-8 




y 


n 


27 


-0.55 


II 


1.1 


8 


K6'!'±l 


4200+150 


0.9 


0.92 


1.67 


3.07E+6 


0.2 


0.18 


l.lE-8 




n 


n 



Notes. Indicates whether the YSO drives a known jet or CO outflow (Morgan et al. 1991; Dent et al. 1998; Allen & Davis 2008; Davis et al. 2009, and this paper). '''' Indicates whether jet lines 
have been detected in our spectra. Used reference/method: from Strom et al. (1989); from Fang et al. (2009); obtained with SpTclass; obtained with the H2O broad band features - water 
breaks; <»' SED modelling. Pre- and outburst values. Fang et al. (2009): L,„i = 1 .9 Lq, SpT = K5, M, = 0.9 Mq, Age = 7E+5 yr, M„,,^ = 1 .5 E-8 Mq yr' . <'''" Fang et al. (2009): Lm = 3.6 Lq, 
SpT = K4.5, M. = 0.9 Mq, Age = 3E+5yr, M„rr=1.9E-8 Mq yr"' . Fischer et al. (2010): Li,„i = 23 L©, M. = 2.2Mo, Age = 3E+5yr, M,,^^ = 4E-7Moyr"'. Strom et al. (1989): SpT=G5-G7. 

Alecianet al. (2009): L. = 97.7 Lq, SpT = B9, M. = 1.9 M©,/?. = 3%, Age = 2E+6yr. Manoj et al. (2006): L. = 93 Lq, SpT = Al, M. >4.9Mo. Hillenbrand et al. (1992): SpT = B9, M. = 
3.6 Mq, R. = 2.8 Rq, M„„ = 5E-6Moyr-'. Fang et al. (2009): Lm = 4.2 L©, SpT = MO, M. = 0.4 M©, Age = 3E+4yr, M„„ = 1.5 E-6M© yr"'. Fang et al. (2009): Lb„i = 60.7 L©, SpT = 
F7, M. = 3Mo, Age = 1.9E+6yr. '■ff> Fang et al. (2009): = 4.9 L©, SpT = Kl, M. = 1.9 M©, Age = 1.4E+6yr, M„„ = 4E-7Moyr-'. 
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5. Discussion 



On the basis of the results from our spectral and photometric 
survey, we can characterise both stellar and cirumstellar proper- 
ties of our flux-limited sample, addressing several issues related 
to the evolution of very young and embedded stellar objects, in 
particular their accretion and ejection activity. 

5.1. Accretion properties of tine sampie 

5.1.1. Accretion luminosities 

Figure 12 shows the Lg^c values plotted as a function of L, for the 
whole sample. The various symbols indicate the different class 
of object, as derived from the standard SED classification (see 
Sect. 4.4.2 and Tab. 7), with Class I, Flat, and Class II repre- 
sented as dots, squares, and triangles, respectively. As already 
observed by several authors (e.g., Muzerolle et al. 1998; Natta 
et al. 2006; Antoniucci et al. 201 1), Lgcc increases with L*. In our 
sample there is a large scatter of Lace I values for objects with 
similar L« values, i. e. about one order of magnitude, excluding 
the outbursting source #25. This spread is similar to what we 
observe in paper-1 for the Cha I and II sample, which have the 
same selection and observational constraints like ours. As a com- 
parison, the Natta et al. (2006) sample of Class II and III YSOs 
in Ophiuchus shows a larger scatter, ranging over more than two 
orders of magnitude. 

Excluding the outbursting source #25, which stands out 
from the others with Lacc/Lt ~ 10, the Lacc/L, ratio for 
our sample ranges from ~0.1 to ~1, indicating that none of 
the observed objects, except source #25 and, marginally, #20 
{L„cc/Lhoi ~0.53), is accretion dominated (i.e. has Lace > L,). 
This result confirms previous findings in other star forming re- 
gions (see e.g. Muzerolle et al. 1998; White & Hillenbrand 
2004; Nisini et al. 2005a; Antoniucci et al. 2008), which indicate 
that several Class I sources are not in the main accretion phase. 
Moreover, there is no clear trend between the Lace I ratio and 
the YSO class, with the various classes randomly distributed be- 
tween 0. 1 and 1 , indicating that high and low accretors are homo- 
geneously present in Class I, Flat, and Class II YSOs. It is worth 
noting that three out of the four sources with Laec/Lt < 0.1 (all 
except source #2) do not show ejection activity in our spectra 
or images (Columns 16 and 17, Tab. 7), indicating that accretion 
and ejection activity are directly related, and confirming the ex- 
pected trend that the higher the accretion rate the higher is the 
ejection rate. 

5.1.2. Mass accretion rates and accretion evolution 

The top panel of Figure 13 shows the mass accretion rates of 
the sampled objects versus their stellar masses. Again, the var- 
ious symbols indicate the different classification of the objects 
(Class I, Flat, and Class II are represented as dots, squares, and 
triangles, respectively), as also reported in Fig. 12. The presence 
of jets/outflows in the images, lines tracing jets in the spectra, or 
both is indicated by open squares, circles, or triangles, respec- 
tively, over-plotted on each object. 

Several authors have pointed out the existence of an empir- 
ical relationship between the mass accretion rate and the stel- 
lar mass, i. e. Maee M^, with k ranging from ~2 (Calvet et al. 
2004; Muzerolle et al. 2003; Natta et al. 2004, 2006) to ~3 (Fang 
et al. 2009). Our data marginally indicate such a trend, because 
of the large spread in Maee (more than 3 orders of magnitude) and 
YSO age, the small M, range, and the small number of objects. 
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thus our data do not allow us to obtain a better estimate for k. 
In the figure we also show the Maee °^ relationship derived 
by Natta et al. (2006) (dashed line) for their Ophiuchus sample, 
and the Maee °^ M;-^ relationship from Fang et al. (2009) (dotted 
line) for their L 1634 and L 1641 samples. The lines serve as a 
guide to the eye and are not intended to fit the data. 

As mentioned in Sect. 4.4.6, the outbursting source #25 
shows the highest mass accretion rate, and, in general, when 
considering a power-law relationship. Class I objects have mass 
accretion rates higher than the other YSOs, being positioned 
mostly in the top-left area of the panel. The remaining Class I 
YSOs are located in the central part of the plot, together with the 
majority of Flat YSOs. Finally, moving towards the bottom-right 
part of the panel, there are Class II and a few Flat YSOs. The bot- 
tom end is occupied by YSOs, that do not show any evidence of 
outflow activity. In conclusion, there is a marginal trend between 
the mass accretion rate evolution and the YSO class. Indeed, this 
is related to the fact that YSO class is not strictly related to YSO 
age, as already shown in Sect. 4.4.5 and Fig. 1 1 

On the other hand, there is a more clear correlation be- 
tween the mass accretion evolution and the age of the YSOs, 
as shown in the bottom panel of Figure 13, which displays the 
same Maee vs M, plot, with different colours indicating different 
age ranges. In this case a clear age stratification is visible in the 
figure, going from the top-left (youngest sources) to the bottom- 
right (oldest sources) of the panel. Moreover, if we consider ob- 
jects in the same age range (i. e. YSOs with the same colour in 
the plot), we note some spread, with high accretors (i. e. with 
large Laee/Lboi ratio) usually positioned at the top-left of their 
group, and low accretors at the bottom-right. For example the 
black points (age<10^ yr) have source #25 with Laee/Lboi ~ 0.9 
at one end, and source #13 with Lgee/Lhoi ~ 0.1 at the other. 
Notably, this low-accreting Class I source does not show any rel- 
evant ejection activity in our images and spectra. The violet sub- 
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Fig. 13. Top panel: Mass accretion rate versus stellar mass for the all 
sample. The different symbols indicate the class of the objects: Class I, 
Flat, or Class II (dots, squares, and triangles, respectively). The object 
ID is also reported. The presence of jets/outflows in the images, lines 
tracing jets in the spectra, or both are indicated by open squares, cir- 
cles, or triangles, respectively. The dashed line shows the Mn„. oc M,' ** 
relationship obtained by Natta et al. (2006) for the Ophiuchus sample, 
whereas the dotted line displays the Mn„. tx M^ ** relationship obtained 
by Fang et al. (2009) for the L 1634 and the L 1641 samples. The lines 
are merely indicative and do not fit the data. Bottom panel: Same as top 
panel, with colours indicating the different age of the objects, ranging 
from <10'^ to 10' yr. 



sample (age -2-5x10 yr) has source # 14 with Lacc/Lboi ~ 0.4 
and source #2 with Lacc/Lh„i ~ 0.1, and so on. This spread is 
usually less than one order of magnitude for the various con- 
sidered age-bins, except for the black points with the outburst- 



ing source #25 (~ 2 orders of magnitude). Although the largest 
spread is due to the outburst, the small spreads (< 1 dex) might 
be attributed to several causes: i) uncertainties on both Mace and 
M, values; ii) uncertainties on the stellar age; Hi) short-time 
scale variability; iv) differences in the initial and environmen- 
tal conditions (see, Natta et al. 2006, for a detailed discussion). 
As stated in Sect. 4.4.6, Mace errors can be up to one dex and 
might cause the spread. There are large uncertainties on the de- 
termination of the stellar age, due to our uncertainties on L, and 
Teff- On the other hand. Tab. 2 indicates that, with the excep- 
tion of the outbursting source #25, short-time scale variability 
cannot be responsible for such a large spread. Indeed, the maxi- 
mum ~ 10 yr variability in the Ks band is <0.5 mag, and, on aver- 
age, <0.2mag. This type of variability is extremely common in 
YSOs (see, e. g.. Beck 2007). Even assuming that extinction has 
not changed and the observed variability in the Kg band is en- 
tirely due to Bry variation (i. e. to accretion), this translates into 
a Bry flux change of -1.6, i.e. less than a factor of 4 in Mace, 
that is well below the desired one dex. 

Despite the large uncertainties on Mace, M,, and age. 
Figure 13 (bottom panel) indicates that the large scatter can 
be, at least partially, explained in terms of YSO age and long- 
time scale variability. Mass accretion evolution has been stud- 
ied for T Tauri discs (Hartmann et al. 1998), and mass accre- 
tion rates are thought to decrease with time as Maee f con- 
sistently with the evolution of a viscous disc (Hartmann et al. 
1998; Muzerolle et al. 2000; Sicilia-Aguilar et al. 2006, 2010). 
According to these viscous disc models, the main infall process 
from envelope to disc is over, thus the envelope is not providing 
a significant amount of mass and angular momentum to the disc. 
Additionally, from this stage on, the quantity of angular momen- 
tum removed by the jet/wind is considered irrelevant. Therefore, 
the disc evolution is mostly driven by viscous processes, which 
control the mass accretion rate. Once the initial disc mass (M^/) 
and viscosity parameter are set, these models predict that the 
mass accretion decreases exponentially with time. Top panel of 
Figure 14 shows the accretion rates of the sampled objects as a 
function of the stellar age, with different colours indicating di- 
fferent mass ranges. Again, the data show a clear trend of the 
mass accretion decaying with time. This trend matches well the 
results of other works (Hartmann et al. 1998; Sicilia-Aguilar 
et al. 2006, 2010), with our data having a smaller Maee spread 
for a given mass interval. For comparison, the black dashed line 
shows the fiducial model of Hartmann et al. (1998) viscous disc. 
The model assumes a stellar mass M» of 0.5 Mq, a disc temper- 
ature at 100 AU of 10 K, an initial disc mass of 0.1 Mq, a 
viscosity parameter a = 10"^, and viscosity exponent j - 1, i. e. 
77 = 1.5 in the Maee f relationship. The curve reproduces 
weU the distribution of our YSOs with 0.4 < M, < 1.2 Mq (i. e. 
violet and red dots). Additionally, more massive and less mas- 
sive YSOs are positioned above and below the curve, respec- 
tively. This is expected, since we observe a coiTelation between 
M„cc and M,. The Hartmann et al. (1998) model would prob- 
ably better match our data using a slightly smaller initial disc 
mass, which would shift the curve downwards, and a slightly 
smaller rj value, which would decrease the curve steepness {-?] 
is the curve slope). This last parameter is particularly important, 
because it regulates the accretion rate, giving us indication of 
the disc lifetime. To get a better estimate of the 77 parameter, we 
then fit those data points with age t > 10^ yr, i.e. in a tempo- 
ral range where the viscous processes likely control the mass 
accretion rate. Additionally, we constrain our sample to stellar 
masses 0.4 < M« < 1.2M0, to be consistent with models and 
data from the literature. The red dotted line shows the best lin- 
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ear fit to the considered data points (namely violet and red dots 
with? > 10^ yr), giving 77 = 1.2+0.4. An identical result (?; ~ 1.2) 
was obtained by Sicilia-Aguilar et al. (2010) from their Cep OB2 
sample. The authors remark that the mass accretion decay ap- 
pears to be slower than previously assumed in the models, where 
77 values range from 1.5 to 2.8 (Hartmann et al. 1998). To take 
into account the Mace dependence from M» and reduce the scat- 
ter in the plot, we normalise Mace dividing by M;. The results 
are shown in the bottom panel of Figure 14. The data scatter is 
smaller, although the most massive stars are still positioned well 
above the average. Considering the data points with f > 10^ yr, 
we obtain 7;= 1.0+0. 3, or, discarding the most massive objects 
(M» >1.5 Mq, i. e. yellow and green data points), 77=1. 3+0.2. 

Finally, on the basis of our previous results, it is reasonable 
to question whether the scatter observed in the Lace vs L* plot of 
Figure 12 might be related to the inferred mass accretion evolu- 
tion. Naively, we might expect an age stratification in the plot 
similar to what we observe in the Mace vs M,, with younger 
objects showing larger Lacc/Lt ratios and the older YSOs with 
smaller ratios. However, if we plot Lace vs L, as a function of 
the age as in Fig. 13, bottom panel, there is not any stratifica- 
tion (see Fig. 15, top panel). Besides the uncertainties on and 
Lace, which indeed affect our results, this is because Lace does 
not depend only on Mace, but also on the M^/Rt ratio, which is 
not linear 

To clarify this point, we simulate the temporal evolution of 
Lace vs L,, for YSOs with 0.4, 0.7, and 1 Mq stellar masses, 
using Siess et al. (2000) evolutionary tracks, and assuming an 
Mace evolution with Mace °^ r^'^- The bottom panel of Figure 15 
shows the results of our simulation, with different colours and 
labels indicating various age bins, from 10* to 10^ yr. According 
to this plot, objects with the same age but with different stellar 
mass usually have different Laee/L, values, which increase with 
Mt. More notably, for a given M,, the Laec/L, ratio does not 
decrease monotonically with time, but, between 10^ and 10^ yr, 
the function shows a local minimum (coloured red in the plot), 
which appears at earlier times with increasing stellar mass (see 
Figure 15, bottom panel). Indeed, our oversimplified simulation 
does not take into account other important factors like mass ac- 
cretion variability, and different environmental and initial con- 
ditions, which may strongly affect the Lgcc value, producing a 
more confused outcome. Nevertheless, it might explain, at least 
partially, the Lacc/L^ spreads so far observed in the literature and 
here. 
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5.1.3. Mitigating the luminosity problem? 

Several observational studies have already shown that a signif- 
icant discrepancy exists between the accretion rates observed 
in embedded YSOs and those predicted by theoretical models. 
This inconsistency was first noted by Kenyon et al. (1990), and 
it is also known as the luminosity problem (Kenyon & Hartmann 
1995), because the early- YSO luminosities, which should be ac- 
cretion dominated, appear to be fainter than expected. The latest 
Spitzer Space Telescope surveys (see, e.g., Evans et al. 2009; 
Enoch et al. 2009) support these findings, revealing that more 
than 50% of the embedded YSOs have bolometric luminosi- 
ties, and thus inferred mass accretion rates, considerably lower 
than those predicted by the free-fall theoretical models (i. e. 
Mace «lQ'^MQyr^ vs ~2xlO"^M0yr-i, for a ClassI YSO 
with M~l Mq; Shu 1977; Terebey et al. 1984). More evidence 
comes from recent NIR spectral studies of Class I YSOs (White 
& Hillenbrand 2004; Nisini et al. 2005a; Doppmann et al. 2005; 
Antoniucci et al. 2008; Beck 2007; Prato et al. 2009), aimed 



Log(0(yr) 



Fig. 14. Top panel: Mass accretion rate as a function of stellar age for 
the observed sample. The various colours indicate different mass values. 
The dashed line displays one of the viscous models of Hartmann et al. 
(1998), for a M, = 0.5 Mq, an initial disc mass = 0.1 Mq, viscosity 
parameter a = 10"-, and viscosity exponent y = 1, i.e. 77 = 1.5 in the 
Mac-c °^ t^'' relationship. The red dotted line shows the best linear fit for 
data points with 0.4 < M. < 1.2 Mq and t > 10^ yr, i.e. M„a. cc r'-^. 
Bottom panel: Same as top panel, but with the M^,^^ values normalised 
to the stellar MJ, to take into account M^cc dependence from M,. 



at measuring both stellar and accretion luminosities, and at es- 
timating the stellar age. These studies have shown that, on av- 
erage, a large fraction of Class I YSOs has mass accretion and 
ejection rates lower than theoretically expected (from one to two 
orders of magnitude), i. e. more similar to the accretion rates of 
the classical T Tauri stars (~10^yr). To reconcile theory with 
observations, these studies provide various motivations. 
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Fig. 15. Top panel: Lace vs L, as in Fig. 12, with colours indicating the 
different age of the objects as bottom panel of Fig. 13. Bottom panel: 
Simulation of the L^,.,. vs L, temporal evolution for YSOs of different 
masses (0.4, 0.7, and 1 Mq), derived considering the evolutionary tracks 
from Siess et al. (2000) and a viscous disc evolution with M^cc r'', 
with 77=1.2. Different colours indicate different age, ranging from 10* 
to 10' yr. 



sample, i. e. three (namely # 6, # 19, and # 20) out of nine Class I 
YSOs have ages ~1 Myr, i. e. they are likely edge-on CTT stars. 

However, results from this and other works indicate that sev- 
eral Class I YSOs are genuine young stellar objects (<5xlO^ yr) 
with accretion rates ranging from ~Macc ~10"^ to 10"** Mq yr"' . 
Therefore, to reach their final masses, these YSOs cannot accrete 
mass in a steady way, but they need to acquire most of their mass 
through short episodic outbursts, during which Mgcc increases by 
some order of magnitudes up to Mfl„.~10"'' Mq yr"' for YSOs 
with ~1 Mq. 

The results of our survey, although biased by the relatively 
small number of targets and by the sample selection criteria, 
are consistent with this picture. A small number of very young 
objects ('true' Class I YSOs) shows high mass accretion rates. 
Among them we detected an outbursting YSO (source #25), 
which has recently increased its mass accretion rate by ~2 orders 
of magnitude (Caratti o Garatti et al. 201 1), whereas the remain- 
ing sources have, on average. Mace ~5xlO"^ Mq yr"' (six ob- 
jects with M, <1 Mq) and Mace ~10"^ Mq yr"' (one object with 
M, ~0.1Mo). Considering different ranges of masses, these 
Maee valucs are almost one order of magnitude higher than those 
reported by White et al. (2007), implying that these young ob- 
jects should spend ~5% of their time in an outburst state to reach 
typical CTT masses, i. .e. this percentage is slightly smaller than 
White et al. (2007) estimates. Therefore, statistically, we should 
observe that among 'true' Class I sources about ~5% of these ob- 
jects are in a outburst phase, while the remaining are quiescent, 
and thus 'under-luminous'. So far our statistics is pretty limited 
and a larger number of deep-NIR imaging surveys is needed. 
Nevertheless, the results of the latest NIR surveys, which dis- 
covered several embedded outbursting YSOs (see, e.g.. Covey 
et al. 2011; Kospal et al. 2011; Caratti o Garatti et al. 2011), 
seem to point in this direction. 



6. Conclusions 

As part of the POISSON project (Protostellar Optical-Infrared 
Spectral Survey On NTT), we present the results of a multi- 
wavelength spectroscopic and photometric survey of 27 embed- 
ded YSOs in the L 1641 star forming region, aimed at deriving 
the stellar parameters and evolutionary stage, as well as inferring 
their accretion and ejection properties. Our multi-wavelength 
database includes low-resolution optical-IR spectra from A^TT 
and Spitzer (0.6-40/^m), as well as photometric data covering a 
spectral range from 0.4 to 1 100/^m, which allow us to construct 
the YSOs spectral energy distributions (SEDs) and to charac- 
terise the object parameters (visual extinctions, spectral types, 
accretion and bolometric luminosities, mass accretion and ejec- 
tion rates). 

The main observational results of this work are the follow- 
ing: 



For example, misclassification might be one of the causes, 
and some Class I objects are likely to be misclassified Class II 
YSOs, seen through optically thick discs. Indeed, geometrical 
effects can easily modify the SED shape and thus the observed 
spectral index (see, e. g., Whitney et al. 2003b). Nevertheless, ra- 
diative transfer models of Whitney et al. (2003b) clearly indicate 
that misclassified objects are about five times fainter than real 
Class I YSOs, and we actually observe these difference in our 
sample (see, e. g., source # 19 vs #23 in Table 7). White et al. 
(2007) suggest that at least one-third of Class I YSOs in Taurus 
have been misclassified, and we find the same result in our small 



- Photometric variability of l^niK >0.3mag is detected in 
-22% of the YSOs (~ 10 yr baseline). 

- The optical-NIR spectra show several emission lines, mostly 
associated with YSO accretion and ejection activity (e. g.. 
Hi, Can, Hei, [Sn], [Feii], H2), as well as other emission 
lines (e.g. Oi, Fez Nai, Mgii, Ci) characteristic of active 
young stars, broad-band molecular absorption features, typ- 
ical of cool objects (VO, TiO and H2O bands) and CO over- 
tone bands both in absorption and emission. Also Spitzer 
spectra show emission lines, mostly H2 pure rotational lines, 
[Feu], [Sin] and [Sim], and [Siii], largely originating from 
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shocks along the YSOs jets, or tracing the disc. Ice absorp- 
tion features at 5 to 8 ;um (H2O) and 15.2yL(m (CO2), as well 
as the amorphous silicate absorption feature at 9.7 /^m have 
been also detected in some of the sources. 

- SEDs were constructed using both photometric and spectro- 
scopic data, allowing to classify the 27 targets as nine Class I, 
eleven Flat, and seven Class II YSOs. The inferred Lhoi val- 
ues range from ~0.3 to ISOL©. 

- The spectral types for the sample range from B7 to M7.5, 
with 74% of the objects being low-mass YSOs (SpT>K5). 
The average and median age of the sample are about 2 and 
1 Myr, respectively, with values ranging from <10^ to 10^ yr 
The average ages of Class I YSOs, Flat, and Class II YSOs 
are 0.4, 2, and 4 Myr Six out of nine Class I YSOs have an 
age of ~10^, while the remaining three appear to be older 
(5x10^ to 10* yr). Also among the Flat sources there are 
three objects (out of eleven) with age >5xlO*yr, i.e. con- 
sistent with more evolved sources. 

- Accretion luminosities, obtained by averaging luminosi- 
ties from Pa^, and Bry lines, range from ~0.1 to ~60Lo' 
whereas mass accretion rates range from 3.6x10"^ to 
1.2x10-5 M0 yr-i. 

From these results we draw the following conclusions: 

- Geometrical effects can significantly modify the SED 
shapes, and the standard SED classification {a computed be- 
tween 2.2 and 24 yum) may not properly reflect YSO age. 

- Mass accretion and ejection rates, extinction, and spectral 
indices decrease with YSO age. The youngest YSOs have 
the highest mass accretion rates, whereas the oldest YSOs 
do not show detectable jet activity in both images and spec- 
tra. Nevertheless, apart from the outbursting source # 25 and, 
marginally, source # 20, none of the remaining YSOs is ac- 
cretion dominated (i. e. having Lace > Lt). 

- There is a marginal trend between the mass accretion rate 
evolution and the YSO class, whereas there is a more clear 
correlation between the mass accretion evolution and the age 
of the YSOs. Despite the large uncertainties on Mace, 
and age, our analysis indicates that the large spread of mass 
accretion rates can be mostly explained in terms of YSO age 
and long-time scale (non-periodic) variability. 

- We obtain a clear correlation between Mace and age (Mace ^ 
r'', with 77=1.2), for YSOs with t > 10^ yr, consistent with 
the mass accretion evolution in viscous disc models. These 
results indicate that the mass accretion decay appears to 
be slower than previously assumed in the classical mod- 
els (Hartmann et al. 1998) (where the rj value ranges from 
1.5 to 2.8) and coincides with the latest results. 

- The majority of very young objects in our sample ('true' 
Class I YSOs, t <5xlQ^ yr) shows, on average, high mass ac- 
cretion rates. Among them, we detected an outbursting YSO 
(source #25), which has recently increased its mass accre- 
tion rate by ~2 orders of magnitude, whereas the remaining 
sources have, on average, Mgcc ~5x 10"^ Mq yr"' (six ob- 
jects with M, <1 Mq) or Mace ~10"^ Mq yr"' when consid- 
ering less massive objects (M» ~0.1 Mq). Considering diffe- 
rent ranges of masses, these Mace values are almost one order 
of magnitude higher than those previously studied, e. g. in 
Taurus and Auriga star forming regions, implying that these 
young objects should spend 1-5% of their lifetime in an out- 
burst state, to reach typical CTTs masses. 
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Table 8. Statistics for the all sample and selected subsamples. The outbursting source #25 has been excluded from any statistic. 



ID 


Avg. (M„,,) 


Avg.(M„,„) 


Avg.(M„„,/M„„) 


Avg. (Age) 


Avg.(L„„/L,„,) 


Avg.(Av) 


Avg. (a) 


Avg.(M.) 


(n)" 


10-' (Moyi-i) 


10-* (M0yi-')(n)° 




(Myr) 




(mag) 




Mo 


All sample (26) 


2.5 


0.29(13) 


0.01 


2.3 


0.26 


12 


0.1 


1.1 


Class 1 (8) 


3.8 


0.34 (5) 


0.01 


0.4 


0.3 


21 


0.7 


0.6 


Flat (11) 


2.6 


0.27 (6) 


0.01 


2.5 


0.25 


11 


0.1 


1.2 


Class 11 (7) 


1.3 


0.20 (2) 


0.03 


4.2 


0.25 


4 


-0.6 


1.6 


Outflow YS0s''(18) 


3.4 


0.29(13) 


0.01 


1.4 


0.28 


13 


0.1 


1.1 


Jetless YSOs' (8) 


0.4 


N/A 


N/A 


4.4 


0.19 


12 





1.0 


Very young YSOs'' (7) 


4.5 


0.32 (6) 


0.01 


0.2 


0.28 


17 


0.7 


0.5 



Notes. Number of elements per subsample This subsample includes the YSOs with an outflow/jet signature in the images (Col. 16 of Tab. 7) 
and/or in the optical/NIR spectra (Col. 17 of Tab. 7); source #25 was not included. This subsample includes the YSOs with no outflow/jet 
signature. This subsample includes the YSOs with age <5xlO' yr (Col. 11 of Tab. 7); source #25 was not included. 
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Appendix A: Photometry & Spectroscopy 



Table A.l. Available photometiy from 0.44 /jm to 8 fim (Jy) 



Source 


B 


V 


R 


Ib 


Ij 


J 


H 


Ks 


IRAC 1 


IRAC 2 


IRAC 3 


IRAC 4 


ID 


0.44/jm 


0.55//m 


0.64//in 


0.79/jni 


0.9//m 


1.25/jm 


1.65//m 


2.2//m 


3.6/im 


4.5 /jm 


5.8//in 


8.0//m 



1 






0.003±0.001 


0.004±0.001 


0.021 ±0.005 


0.0224±0.0009 


0.066±0.002 


0.089±0.002 


0. 122±0.007 


0.1 15±0.006 


0.103±0.006 


0.097±0.005 


2 


6e-5±2e-5 




0.001±0.0003 


0.005±0.001 


0.0055±0.0005 


0.048±0.001 


0.166±0.005 


0.347±0.007 


0.651±0.006 


0.83±0.001 


0.97±0.01 


1.17±0.01 


3 


0.12±0.03 


0.23±0.05 




0.43±0.06 


0.25 ±0.03 


0.91 ±0.02 


1.68±0.04 


2.79±0.06 


2.09±0.09 




4.1 ±0.1 


5.0±0. 1 


4 


0.0016±0.0003 


0.0043±0.0008 


0.007 ±0.0019 


0.018±0.003 


0.018±0.003 


0.059±0.002 


0.089±0.003 


0. 104±0.003 


0. 116±0.002 


0.137±0.004 


0.156±0.004 


0.221±0.006 


5 






4e-5±le-5 


23e-5±6e-5 


0.0004±0.0002 


0.0092±0.0003 


0.03911 ±0.0009 


0.076±0.002 


0.147±0.003 


0.203±0.004 


0.245±0.007 


0.382±0.007 










<9e-5 


2e-5±le-5 


0.00 10±0. 0005 


0.004±0.002 


0.0128±0.0008 


0.097±0.004 


0.184±0.003 


0.240±0.007 


0.266±0.005 


7 


0.003±0.002 




0.028±0.001 


0.045 ±0.002 


0.19±0.005 


0. 169±0.003 


0.265±0.007 


0.43±0.01 


0.73±0.04 


0.86±0.05 


0.88±0.05 


1.70±0.09 


g 


3e-5±le-5 




0.0006±0.0002 




0.00024±0. 00007 


0.317±0.007 


2.02±0.06 


5.8±0. 1 






15.7±0. 1 




9 






0.010±0.001 


0.022±0.002 


0.032±0.005 


0. 128±0.003 


0.289±0.007 


0.41±0.01 


0.372±0.007 


0.324±0.006 


0.274±0.007 


0.323 ±0.006 


10 






0.00023±0. 00004 


0.00095±0. 00005 


0.0008±0.0003 


0.0 123±0. 0003 


0.0290±0.0006 


0.040±0.001 


0.084±0.005 


0.108 ±0.008 


0.14±0.01 


0.19±0.01 


11 






0.005 ±0.002 


0.013±0.002 


0.050±O.O03 


0. 105±0.002 


0.220±0.003 


0.36±0.01 


0.44±0.03 


0.49±0.03 


0.55±0.03 


0.52±0.03 


12 






35e-5±5e-5 


0.0010±0.0002 


0.0004±0.0002 


0.0123±0.0006 


0.0220±0.0006 


0.0242±0.0006 


0.021 ±0.001 


0.0236±0.001 


0.032±0.002 


0.050±0.003 


13 












0.00029±0.0001 


0.0007 ±0.0002 


0.0027±0.0009 


0.054±0.005 


0.08 1± 0.006 


0.075±0.006 


0.090±0.007 


14 












0.0010±0.0001 


0.0089±0.0002 


0.042±0.001 


0.241 ±0.004 


0.413±0.004 


0.576±0.009 


0.793±0.009 


15 


9e-5±2e-5 




0.0009±0.0002 


0.00 18±0. 0005 


0.0031 ±0.0005 


0.0093±0.0002 


0.0285 ±0.0007 


0.055±0.001 


0.099±0.003 


0. 116±0.003 


0.127 ±0.005 


0.177±0.003 


16 


0.0013±0.0003 




0.0156±0.004 


0.04±0.01 


0.045±0.005 


0.202±0.004 


0.329±0.008 


0.45±0.02 


O.507±O.O05 


0.490±0.009 


0.453±0.008 


0.774±0.007 


17 












37e-5±4e-5 


0.0084±0.0002 


0.0410±0.0009 


0.102±0.002 


0.163±0.004 


0.230±0.006 


0.325±0.006 


18 


0.012±0.004 


0.022±0.008 






0.07±0.01 


0.091 ±0.002 


0.120±0.0O3 


0.129±0.003 


0.19±0.01 


0.26±0.01 


0.34±0.02 


0.52±0.03 


19 










2e-5±le-5 


5e-5±le-5 


0.0004±2e-5 


0.0076±0.0002 


0.036±0.001 


0.049±0.001 


0.060±0.003 


0.080±0.003 


20 












0.0003±0.0001 


0.0060±0.0002 


0.0399±0.0007 


0.218±0.004 


0.314±0.006 


0.419±0.008 


0.474±0.008 


21 












<0.0003 


0.0013±0.0001 


0.0121±0.0003 


0.043± 0.002 


0.069±0.001 


0.087±0.003 


0.099±0.003 


22 










2e-5±le-5 


0.0013±6e-5 


0.0079±0.0002 


0.0158±0.0004 


0.0106±0.001 


0.020±0.002 


0.022±0.002 


0.025±0.006 


23 










2e-5±le-5 


0.0008±6e-5 


0.0087±0.0002 


0.0386±0.0004 


0.136±0.003 


0.253±0.005 


0.352±0.007 


0.45±0.01 


24 










3e-5±le-5 


86e-5±6e-5 


0.0056±0.0002 


0.0250±0.0006 


0.090±0.002 


0.167±0.003 


0.302±0.008 


0.386±0.007 


25 










le-5±le-5 


0.00044±4e-5 


0.0040±0.0001 


0.0132±0.0004 


0.121 ±0.007 


0.18±0.01 


0.20±0.01 


0.23±0.01 


26 








9e-5±le-5 


0.0002±le-5 


0.0107±0.0002 


0.061±0.001 


0.154±0.003 


0.22±0.01 


0.29±0.02 


0.32±0.02 


0.32±0.02 


27 


5e-5±le-5 




30e-5±7e-5 


0.001 1±0.0003 


0.0019±0.0009 


0.0120±0.0002 


0.0324±0.0006 


0.0479±0.0009 


0.042±0.002 


0.045±0.002 


0.030±0.002 


0.024±0.002 



Table A.2. Available photometry from 24 fim to 1 100 /im (Jy) 



Source MIPS 1 MIPS 2 MIPS 3 UKT14 SCUBA UKT14 SCUBA UKT14 
ID 24 /jm IQiim 160/im 350;jm 450;jm 800 ^m 850 /jm llOO/im 



1 


0.172±0.005 












2 


4.5±0.3 


11.2±0.5 




6.52±2.1 


3.75±0.27 


0.635±0.032 


3 


6.6±0.8 


6.5±2.2 






2.540±0.06 


0.9±0.06 


4 


0.88±0.08 


1.5±0.2 










5 


0.97±0.04 


1.4±0.1 


2.5±1 








6 


1.2±0.1 


5.2±0.6 


4±1 




1.42±0.13 


0.19±0.014 


7 


2.6±0.1 


4.1±0.2 


6±1 








8 


52±2 


97±1 


35±9 


13.7±0.29 


9.04±0.06 


1.90±0.087 


9 


0.74±0.02 


0.8±0.1 










10 


0.48±0.05 


0.6±0.1 










11 


1.32±0.07 


4.0±0.7 


<6 








12 


0.45±0.03 


0.8±0.2 










13 


1.7±0.1 


5.4±0.9 


5.3±2.2 








14 


2.3±0.2 


5.4±0.2 


5±1 








15 


O.51±0.02 


1.6±0.2 


<4 






<0.2 


16 


3.2±0.2 


1.6±0.2 


<0.8 






<0.2 


17 


0.74±0.05 


1.3±0.1 


3.5±1 






0.450±0.022 


18 


1.91±0.07 


2.3±0.3 










19 


0.31±0.05 


1.9±0.5 










20 


1.4±0.1 


4.9±0.5 


3±1 






0.148±0.06 


21 


0.26±0.01 


4.4±0.4 










22 


0.1±0.02 












23 


4.5±0.2 


10.7±1.2 


3.5±1.3 




2.3±0.5 


0.7±0.3 


24 


2.9±0.2 


10.7±0.6 


5.5±1.4 




2.79±0.64 




25 


0.73±0.08 


1.4±0.2 


0.5±0.1 




0.2±0.1 




26 


0.26±0.01 


0.5 ±0.2 










27 


0.148±0.008 


0.96±0.22 











0.239±0.039 



1.40±0.06 0.767±0.06 



0.7±0.011 0.289±0.030 



Caratti o Garatti et al.: A muld-wavelength spectroscopic and photometric survey of young protostars in L 1641, Online Material p 3 
Table A.3. Observed emission lines in [CSH2001]13811 - ID 1. 



Line id. 


A 


F 


AF 


FWHM 


W,i 




(yum) 


(erg cm 




A 


A 






EFOSC & Sofl - NTT 






Ua 


0.656 


9.53e-14 


l.OOe-14 


11.0 


-11.01 


Pa 5 


1.005 


8e-15 


3e-15 


16.8 


-2.0 


Ci 


1.069 


7.42e-15 


2.00e-15 


18.0 


-18.0 


Pay 


1.094 


1.09e-14 


0.30e-14 


19.1 


-2.8 


PajS 


1.282 


1.8e-14 


0.30e-14 


21.6 


-4.0 


Bry 


2.166 


2.96e-14 


0.50e-14 


55.0 


-5.1 



Line id. 


A 


F 


AF 


FWHM 


W,i 




(yum) 


(Jyjum) 


(jum) 


(jum) 






IRS 


- Spitzer 






H2 


17.0 


0.010 


0.002 


0.25 


-0.07 


[Fell] 


22.9 


0.014 


0.003 


0.24 


-0.08 


[Fell] 


24.5 


0.007 


0.003 


0.30 


-0.04 


[Sim] 


25.6 


0.028 


0.005 


0.40 


-0.16 


H2 


28.2 


0.031 


0.005 


0.40 


-0.16 


[Sm] 


33.5 


0.053 


0.008 


0.6 


-0.27 


[Fell] 


35.8 


0.023 


0.007 


0.24 


-0.09 


[Fei] 


36.5 


0.068 


0.008 


0.36 


-0.27 
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Table A.4. Observed emission lines in [CTF93] 50 - ID 2. 



Line id. 


A 


F 


AF 


FWHM 








(erg cm ^ s ') 


A 


A 




EF08C & 8ofI - 


NTT 






[Ol] 'D2-3P2 


0.630 


7.66e-15 


0.16e-15 


8.322 


-7.065 


Fell 


0.635 


0.96e-15 


0.30e-15 


11.8 


-0.87 


[Oi] 'D2-^Pi 


0.636 


2.66e-15 


0.30e-15 


9.8 


-2.343 


Fe II 


0.646 


0.78e-15 


0.30e-15 


8.90 


-0.564 


Fe I 


0.649 


1.06e-15 


0.40e-15 


12.80 


-0.857 


Fe II 


0.651 


l.lOe-15 


0.40e-15 


11.00 


-0.6 


[Ni] 'D2-^Po 


0.653 


4.43e-14 


0.50e-14 


11.8 


-1.246 




0.656 


9.87e-15 


0.40e-15 


6.76 


-6.96 


Ha PCygni 


0.655 


-1.6e-15 




4.5 


1.1 


[S n] 205/2-^83/2 


0.672 


3.30e-15 


0.30e-15 


12.08 


-2.14 


[Sn] 203/2-^83/2 


0.673 


5.49e-15 


0.25e-15 


9.71 


-3.4 


He I 3Po2 -0-^81 


0.706 


8.36e-15 


0.30e-14 


36.2 


-4.845 


Fe II 


0.707 


l.Ole-14 


0.50e-14 


12.2 


-0.296 


? 


0.710 


l.OOe-14 


0.30e-15 


35.13 


-6.16 


? 


0.713 


1.27e-13 


l.OOe-14 


12.4 


-4.45 


? 


0.715 


0.93e-14 


0.30e-15 


23.64 


-5.12 


7 


0.733 


1 £.r\^ 1 c 

1.60e-15 


0.4e-15 


14.9 


-0.87 


7 


0.12)1 


2.70e-15 


0.4e-15 


14.8 


-1.424 


Pa 20 


0.839 


2.22e-15 


O.lOe-14 


15.1 


-0.77 


Fa iv 


0.842 


2.48e-14 


l.OOe-14 


111 
11.1 


-2.619 


Do 1 C 1 


r\ A A 

U.844 










Oi^'Po- 8j 
Oi^P2-^8j 


r\ A c 

0.845 


5.55e-15 


O.lOe-14 


11.3 


-2.07 












Oi^Pi-^8* 












Pa 16 + 


0.850 










Can2p3/2-2D3/2 


0.850 


1.86e-14 


0.20e-14 


9.9 


-6.66 


Cail2p3/2-2D5/2 + 


0.854 


1.82e-14 


0.20e-14 


9.9 


-6.42 


Pa 15 


0.855 










Pa 14 


0.860 


2.86e-15 


O.lOe-14 


12.40 


-1.0 




0.860 


2.00e-15 


O.lOe-14 


9.90 


-0.70 


Ni 


0.863 


3.88e-15 


0.20e-14 


20.90 


-1.34 


Cail2Pi/2-2D3/2 


0.866 


1.65e-14 


0.20e-14 


10.15 


-5.70 


Pa 13 


0.867 


4.43e-15 


0.20e-14 


16.41 


-1.52 


Pa 12 


0.875 


4.92e-15 


0.20e-14 


17.9 


-1.70 


Mgi 


0.881 


4.07e-15 


0.20e-14 


18.6 


-1.43 


Pa 5 


1.005 


2.92e-13 


O.lOe-13 


21.78 


-13.1 


Ci(l) 


1.069 


1 *7 1 „ ^ A 

1.71e-14 


0.40e-14 


20.2 


-2.92 


He iPCygni 


i .Uo J 


'7 C^a 1 /I 

-z.oe-i4 




15.8 


4.64 


8i lor Fe 11 


1.087 


5.31e-14 


l.OOe-14 


20.4 


-2.418 


Pay 


1.094 


2.02e-14 


0.40e-14 




i.ll 1 


C I (/4j 


1.175 


1.23e-14 


0.80e-14 


1 "2 


1 


[r nj i'2- U2 or I 


1.189 


2.11e-14 


0.80e-14 


io.U 


/.JO 


8ii? 


1.205 


2.46e-14 


0.80e-14 


21.1 


-2.86 


[Fell] a4D7/2-a«D;^2 
He I? 


1.257 


1.76e-14 


0.70e-14 


26.9 


-1.76 


1.270 


3.61e-14 


0.80e-14 


16.5 


-3.59 


Pa/? 


1.282 


4.5e-14 


0.40e-14 


18.04 


-4.4 


[Fe 11] a'^D7/2-a'^F9/2 


1.644 


2.11e-14 


0.80e-14 


44.0 


-1.156 


Bry 


2.166 


7.41e-14 


0.80e-14 


45.62 


-3.33 



Line id. ^ F AF FWHM 
(jum) (Jyjum) (yum) (yum) 



IR8 - 8pitzer 



[8 III] 


33.5 


0.18 


0.06 


0.4 


-0.24 


[8iii] 


35.3 


0.13 


0.06 


0.2 


-0.17 


[Fell] 


35.8 


0.19 


0.07 


0.3 


-0.24 
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Table A.4. continued. 



Line id. 


A 


F AF 


FWHM 








(Jyjum) 


(jum) 


(jum) 


[Fei] 


36.5 


0.22 0.07 


0.3 


-0.26 
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Table A.5. Observed emission lines in [CTF93] 47 - ID 3. 



Line id. 


A 


F 


AF 


FWHM 








(erg cm 


-2 s-i) 


A 


A 




EFOSC & Sofl - NTT 






Fell 


0.615 


3.17e-13 


O.lOe-13 


9.77 


-1.755 


Fe II 


0.625 


6.69e-13 


O.lOe-13 


14.4 


-3.68 


Fe II 


0.635 


3.13e-13 


O.lOe-13 


13.10 


-1.781 


Fe II 


0.637 


3.17e-13 


O.lOe-13 


13.00 


-1.754 


Fe II 


0.642 


2.32e-13 


O.lOe-13 


11.10 


-1.27 


Fei 


0.643 


3.04e-13 


O.lOe-13 


11.80 


-1.67 


Fe II 


0.646 


6.78e-13 


O.lOe-13 


10.00 


-3.725 


Fei 


0.649 


1.81e-13 


O.lOe-13 


16.6 


-1.585 


Fei 


0.651 


4.56e-13 


O.lOe-13 


10.58 


-2.502 


Hot 


0.656 


1.33e-ll 


l.OOe-14 


9.68 


-73.17 


Hei^Poi -0-^Si 


0.667 


9.35e-14 


l.OOe-14 


11.48 


-0.508 


Fei? 


0.697 


1.26e-13 


2.00e-14 


16.3 


-0.68 


Fe II 


0.746 


2.89e-13 


3.00e-14 


12.8 


-1.76 


Fe II 


0.751 


1.98e-13 


3.00e-14 


12.4 


-1.22 


Fe II 


0.771 


3.64e-13 


4.00e-14 


12.0 


-2.42 


Oi^Po-^s; 


0.777 


5.36e-13 


4.00e-14 


10.79 


-3.567 


Mgii 


0.787 


2.22e-13 


5.00e-14 


17.7 


-1.54 


Mgii 


0.789 


2.72e-13 


5.00e-14 


12.6 


-1.91 


Si I 


0. /y4 


1. /se-13 


C AA^ 1 A 

5.00e-14 


10.6 


-1.285 


Pa 22 


0.836 


1 ^ 1 ") 
1.53e-13 


C AA^ 1 A 

5.00e-14 


14.1 


-1.198 


Pa 19 


0.842 


1 T 1 1 

1.30e-13 


C AA^ 1 A 

5.00e-14 


10.6 


-1.0 


ra io + 


U.o44 










Oi^Pq- Sj 


0.545 


s.//e-13 


C AA^ 1 A 

5.00e-14 


1 1 AA 

1 1.00 


-6.37 


Oi^P2--'Sj 












Oi^Pi-^S* 












Pa 17 


0.847 


1.89e-13 


5.00e-14 


10.0 


-1.46 


Pa 16 + 


0.850 










Cail2p3/2-2D3/2 


0.850 


3.13e-12 


5.00e-14 


10.0 


-24.2 


Caii%2-^D5/2 + 


0.854 


3.44e-12 


5.00e-14 


10.3 


-26.6 


Pa 15 


0.855 










Pa 14 


0.860 


2.80e-13 


0.40e-13 


12.0 


-2.15 


Cail2Pi/2--D3/2 + 


0.866 


3.04e-12 


5.00e-14 


10.32 


-23.67 


Pa 13 


0.867 










Pa 12 


0.875 


3.39e-13 


5.00e-14 


12.7 


-2.77 


Mgi 


0.881 


3.57e-13 


5.00e-14 


14.7 


-3.13 


Fe I 


0.882 


2.20e-13 


5.00e-14 


16.7 


-1.95 


Pa 11 


0.887 


3.63e-13 


5.00e-14 


10.0 


-3.24 


Can2F°^2- D3/2 
Pa 10 


0.891 


2.55e-13 


5.00e-14 


13.0 


-2.3 


0.902 


4.02e-13 


l.OOe-13 


12.9 


-4.7 


Pa 9 


0.923 


8.79e-13 


l.OOe-13 


26.0 


-10.42 


Pa 8 


0.955 


1.06e-12 


0.80e-13 


17.91 


-5.67 


Fe II 


1.000 


1.25e-12 


0.80e-13 


17.0 


-5.35 


Pa 5 


1.005 


1.52e-12 


0.80e-13 


16.00 


-8.24 


Fe II 


1.051 


7.83e-13 


0.80e-13 


20.8 


-5.35 


Ci(l) 


1.069 


5.91e-13 


0.80e-13 


20.50 


-3.19 


Ci(l) 


1.074 


1.33e-13 


8.00e-14 


21.8 


-0.71 


He I 




Q A 1 1 


A Qf\a 1 T 


18.4 


-5.08 


Si lor Fe 11 


1.087 


5.93e-13 


8.00e-14 


17.8 


-3.171 


Pay 


1.094 


2.83e-12 


0.80e-13 


22.7 


-16.2 


[Pii]? 


1.114 


4.91e-13 


8.00e-14 


18.8 


-2.68 


Can P3/2- Si/2 


1.184 


7.63e-13 


1.50e-13 


io.o 


-j./y 


He I? 


1.270 


2.61e-13 


0.80e-13 


19.2 


-1.713 


Pa/? 


1.282 


3.86e-12 


0.80e-13 


22.1 


-23.2 


Mgi3pO-3Si(m) 


1.505 


6.55e-13 


l.OOe-13 


29.4 


-3.18 


Br 16 


1.556 


2.27e-13 


0.80e-13 


38.6 


-1.11 


Nai2Si/2-2p"(m) Mgi 


1.576 


5.53e-13 


0.80e-13 


52.6 


-2.72 


Br 14 


1.588 


6.47e-13 


0.80e-13 


47.00 


-3.12 
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Table A.5. continued. 



Line id. 


A 


F 


AF 


FWHM 






(//m) 


(erg cm 




A 


A 


Br 13 


1.611 


6.29e-13 


0.80e-13 


31.17 


-3.04 


Br 12 


1.641 


9.12e-13 


0.80e-13 


47.35 


-6.905 


Ur 1 1 


i .UO i 




W. Ov7C i J 




-4 Q'? 


Nai2D3/2-5/2-^P°(m)Ci? 


1.689 


6.76e-13 


0.80e-13 


40.8 


-3.3 


7 


1.734 


4.07e-13 


0.80e-13 


29.5 


-2.01 


Br 10 + ? 


1.737 


7.67e-13 


0.80e-13 


31.19 


-3.683 


? 


1.742 


3.58e-13 


0.80e-13 


29.5 


-1.77 


Bry 


2.166 


1.81e-12 


0.80e-13 


35.00 


-13.0 


Nai2D3/2-5/2-'P°(m) 


2.207 


2.13e-13 


0.80e-13 


40.60 


-1.24 
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Table A.6. Observed emission lines in [CTF93] 32 - ID 4. 



Line id. 


A 


F 


AF 




FWHM 






(yum) 


(erg cm 


-2 s-i) 




A 


A 




EFOSC & Sofl - 


IN 1 1 








[Oi] D2- P2 




1 QAp izt 


A 1 

0.15e- 


1 A 

■14 


y./a 




[Oi] 'D2- Pi 


U.DJO 


j.oje-i J 


1.50e- 


■15 


13.90 


-1.164 


Ha 


U.OjO 


o./je-i J 


0.20e- 


■14 


10.4 


-143.1 


He I P02 — U— 01 


U.66 / 


T 1 1 C 

/.i5e-i5 


A OA^ 

O./Oe- 


1 A 

■14 


1 1 CA 

1 1.50 


-1.35 


Do 1 1 

ra is + 


U.o44 












Oi^Po- Sj 


0.845 


1.28e-14 


A A A^ 

0.40e- 


1 A 

■14 


1 A /I 

10.4 


A 
-Z.Z/4 


Ol^P2- Sj 














Oi^Pi-^Sj 














Pa 16 + 


i\ OCA 

O.o50 












Can2p3/2-2D3/2 


0.850 


2.26e-14 


0.40e- 


■14 


11.0 


-4.0 


Cail2p3/2-2D5/2 + 


0.854 


2.26e-14 


0.40e- 


■14 


10.5 


-4.0 


Pais 


0.855 












Caii^Pi/2-^D3/2 + 


0.866 


3.1e-14 


0.40e- 


■14 


11.9 


-5.5 


Pa 13 


0.867 












Pad 


1.005 


l.OOe-14 


0.30e- 


■14 


24.65 


-5.69 


He I 


1.083 


1.53e-14 


0.30e- 


■14 


20.1 


-8.55 


Pay 


1.094 


1.18e-14 


0.30e- 


■14 


20.9 


-8.14 


Paj8 


1.282 


1.65e-14 


0.30e- 


■14 


25.18 


-9.2 


Bry 


2.166 


3.96e-15 


0.13e- 


■15 


41.72 


-4.54 



Line id. 


A 


F AF 


FWHM 






(yum) 


(Jyjum) 


(jum) 


(jum) 






IRS - Spitzer 






[Fell] 


22.9 


0.016 0.006 


0.30 


-0.02 


[Fell] 


24.5 


0.038 0.006 


0.30 


-0.02 


H2 


28.2 


0.025 0.006 


0.38 


-0.02 


[Sm] 


33.5 


0.053 0.01 


0.5 


-0.14 
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Table A.7. Observed emission lines in [CTF93] 72 - IDS. 



Line id. 


A 


F 


AF 


FWHM 


w. 






(erg 




A 


A 



EFOSC & Sofl - NTT 



Ha 


0.656 


1.33e-15 


0.20e- 


■15 


7.75 




Pa 20 


0.839 


1.58e-15 


0.30e- 


■15 


11.9 


-24.5 


Pa 18 


0.844 


1.46e-15 


0.30e- 


■15 


15.14 


-26.3 


Pa 16 


0.850 


4.53e-16 


0.30e- 


■15 


10.7 


-5.7 


Cail2p3/2-2D3/2 


0.850 


4.15e-15 


0.30e- 


■15 


10.6 


-19.0 


Can2p3/2-2D5/2 + 


0.854 


4.41e-15 


0.30e- 


■15 


11.1 


-54.5 


Pa 15 


0.855 


3.77e-16 


0.30e- 


16 


12.1 


-4.5 


Cail2Pi/2-2D3/2 


0.866 


3.84e-15 


0.30e- 


■15 


10.42 


-17.5 


Pa 12 


0.875 


1.51e-15 


0.40e- 


■15 


9.7 


-15.07 


[S n] 2Pi/2-2D3/2 


1.033 


3.08e-15 


0.80e- 


■15 


18.1 


-4.8 


Ci(l) 


1.069 


4.31e-15 


l.OOe- 


■15 


19.30 


-8.16 


He 1+ rCygm 


1.083 


3.28e-15 


0.80e- 


■15 


A 1 


-3.1 


Si lor Fe ii 


1.087 


z.yze-15 


0.80e- 


1 C 

■15 


A 1 

20.1 


-3.2 


Fay 


1.094 


C\'~) ^ 1 c 

8.93e-15 


A A^ 


1 C 

■15 


"^A C 

20.5 


-9.46 


C I (z4) 
Can P3/2- Si/2 


1.1 /j 


4. /Ue-i J 


l.UUe- 


1 c 

■1 J 


1 /I 
iy.4 


-J. 43 


1 ^ A 
1.1 o4 


/I C /I 1 c 

4.54e-i5 


1 AA^ 

l.UUe- 


1 c 

■1 J 


1 T A 
1 /.O 


-3.3 / 


[Pii] "P2- D2 orCi 


1 1 or\ 

1.189 


4.60e-15 


1 AA^ 

l.OOe- 


■15 


22.6 


1 /I 
-3.14 


Mgi 'Pi-'D2 or Si I 


1.198 


3.06e-15 


l.OOe- 


■15 


23.3 


-3.6 


Si I? 


1.205 


3.05e-15 


l.OOe- 


■15 


20.0 


-1.95 


Si I? 


1.218 


1.06e-14 


l.OOe- 


■15 


17.0 


-1.95 


[Feii]a4D7/2-a^D;^2 
He I.' 


1.257 


8.82e-15 


l.OOe- 


■15 


21.7 


-4.6 


1 OTA 


4.0 /e-i5 


i.OOe- 


■15 


1 1 
16.1 


-/.3a 


PayS 


1.282 


2.79e-14 


O.lOe- 


■14 


20.09 


-12.7 


[Fell] ? 


1.321 


6.03e-15 


1.50e- 


■15 


19.0 


-2.53 


[Fen] ? 


1.325 


6.01e-15 


1.50e- 


■15 


24.9 


-2.53 


Mgi¥"-3Si(m) 


1.505 


2.20e-14 


0.40e- 


■14 


25.3 


-5.51 


[Fen] a'*D7/2-a''F9/2 


1.644 


2.90e-14 


2.00e- 


■15 


26.2 


-1.86 


H2 


1.735 


1.16e-14 


0.40e- 


■14 


28.5 


-2.38 


H2 


1.748 


1.63e-14 


0.50e- 


■14 


33.5 


-3.78 


H2 


2.121 


1.40e-14 


0.40e- 


■14 


40.23 


-2.8 


Bry 


2.166 


3.754e-14 


0.40e- 


■14 


43.5 


-7.6 


Nai2D3/2-5/2-'P°(m) 


2.207 


9.00e-15 


0.40e- 


■14 


59.00 


-1.809 



Line id. 


A 


F 


AF 


FWHM 






(fim) 


(Jyyum) 






(p.m) 






IRS - Spitzer 








H2 


28.2 


0.028 


0.009 


0.39 


-0.02 


[Sm] 


33.5 


0.042 


0.01 


0.4 


-0.03 


[Sin] 


34.8 


0.062 


O.OI 


0.4 


-0.04 


[Fe n] 


35.8 


0.050 


0.01 


0.3 


-0.03 
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Table A.8. Observed emission lines in [CTF93] 83 - ID 6. 



Line id. 


A 


F AF 


FWHM 


w,, 




(yL(m) 


(erg cmr'^ s"') 


A 


A 



EFOSC & Sofl - NTT 
H2 1-08(1)" 2.122 2e-15 le-15 36.00 -5.4 
Bry 2.166 4e-15 le-15 82.35 -6.5 



Line id. A ¥ tsP FWHM W,i 
(/^m) (Jy/um) (jum) {^lm) 



IRS - Spitzer 



[Sim] 


25.6 


0.067 


0.02 


0.53 


-0.04 


H2 


28.2 


0.041 


0.01 


0.36 


-0.02 


[SlIl]H- ? 


33.5 


0.12 


0.02 


0.56 


-0.06 


[Sill] 


35.3 


0.14 


0.03 


0.47 


-0.06 
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Table A.9. Observed emission lines in [CTF93] 62 - ID 7. 



Line id. 


A 


F AF 


FWHM 






(yum) 


(erg cmT^ s"') 


A 


A 






EFOSC & SofI - NTT 






Ha 


0.656 


1.18e-13 l.OOe-14 


9.2 


-5.7 


Bry 


2.166 


8.47e-14 1.70e-14 


53.4 


-3.1 



Line id. 


A 


F 


AF 


FWHM 






{^lm) 


(Jy/zm) 


(jum) 


(jum) 






IRS ■ 


■ Spitzer 






[Sim] 


25.6 


0.057 


0.009 


0.36 


-0.02 


H2 


28.2 


0.05 


0.01 


0.37 


-0.02 


[Sm] 


33.5 


0.05 


0.02 


0.32 


-0.06 


[Sin] 


34.8 


0.06 


0.02 


0.4 


-0.03 


[Fei] 


36.5 


0.11 


0.03 


0.5 


-0.04 
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Table A.IO. Observed emission lines in [CTF93] 79 - ID 8. 



Line id. 


A 


F AF 


FWHM 


W,i 




(yum) 


(erg cmT^ s"') 


A 


A 






EFOSC & SofI - NTT 






He I 


2.058 


1.04e-12 0.30e-13 


34.60 


-3.346 


Br 7 


2.166 


1.38e-12 l.OOe-13 


45.42 


-3.639 
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Table A.ll. Observed emission lines in [CTF93] 99 - ID 9. 



Line id. 


A 


F 


AF 


FWHM 






(yum) 


(erg cm" 




A 


A 






EFOSC & SofI - NTT 






Ha 


0.656 


1.34e-13 


0.02e-13 


9.02 


-16.99 


Ci(l) 


1.069 


3.92e-14 


l.OOe-14 


33.2 


-2.097 


He I 


1.083 


3.81e-14 


0.80e-14 


14.95 


-1.93 


He iPCygni 


1.082 


-2.30e-14 




10.2 


1.2 


Ci? 


1.167 


3.11e-14 


l.OOe-14 


19.05 


-1.359 


C I (24) 


1.175 


3.08e-14 


l.OOe-14 


16.0 


-1.349 


7 


1.263 


5.06e-14 


1.20e-14 


36.2 


-1.99 


He I? 


1.276 


3.14e-14 


l.OOe-14 


17.5 


-1.22 


Pay? 


1.282 


4.20e-13 


1.20e-14 


28.03 


-1.64 


He I 


2.058 


5.42e-14 


l.OOe-14 


28.86 


-2.02 


H2 1-0 S(l) 


2.122 


3.05e-14 


l.OOe-14 


28.86 


-1.3 


Bry 


2.166 


9.80e-14 


1.20e-14 


64.21 


-3.79 




Line id. 


A 


F 


AF 


FWHM 








(Jyjum) 


(jum) 


(jum) 






IRS - Spitzer 






[Fell] 


22.9 


0.05 


0.01 


0.37 


-0.04 
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Table A.12. Observed emission lines in [CTF93] 87 - ID 10. 



Line id. 


A 


F 


AF 


FWHM 








(erg cmT^ 




A 


A 



EFOSC & SofI - NTT 







9 A/Lf> 1 S 

z.o^e-i J 


w.zwe-i J 




1 7 9S 
- i / .ZJ 










4. Q 


7 7S 


[UIJ U2 ri 


\J.\JJ\J 


/ .uie-io 


z.uwe-io 


A AC\ 


A 1 4 




635 






3 6 


1 75 


Ha 


0.656 


3.05e-14 


0.03e-14 


10.12 


-96.09 


Her ^'Poo - 


0.667 


6 OOe-16 


3 OOe-16 


9.6 


-2.18 


uj U5/2 03/2 + 


w.u / z 


7?p 1 S 


^Op 1 S 


1 no 


Z.JO i 


LSiij U3/2- 03/2 


U.O / J 


u.ooe-i J 


"^nf. 1 ^ 
u.^ue-i J 


y.ou 


1 9A 


rie I r O2 — U— 01 + 


U. /Uu 


1 AQ« 1 /I 


M 1 /l^:. 1/1 

U. i4e-i4 


1 1 

i i . J 


1 C^9 
-i.OJZ 


Pn 1 S -1- 

-Ta 1 T 


n Sztzt 










T 3t> 3c* 
Ul Vq- 


U.o4j 


1 ^ 
/. jye-i J 


u.oue-i J 


1 n 1 9 
iU. iz 


-J.JJ 


Oi P2- Sj 












Oi^Pi-^SJ 












Pa 16 + 


0.850 










CaiI-P3/2-^D3/2 


0.850 


1.28e-14 


0.80e-15 


10.2 


-18.92 


Can2p3/2-2D5/2 + 


0.854 


1.29e-14 


0.80e-15 


10.4 


-18.64 


Pa 15 


0.855 










Pa 14 


0.860 


1.60e-15 


0.80e-15 


11.80 


-2.237 


Ni 


0.863 


2.70e-15 


0.80e-15 


18.18 


-3.76 


Cail2Pi/2-2D3/2 


0.866 


1.06e-14 


0.80e-15 


9.9 


-14.12 


Pa 13 


0.867 


2.00e-15 


0.80e-15 


17.0 


-2.63 


Pa (5 


1.005 


5.25e-15 


0.80e-15 


15.45 


-3.3 


Ci(l) 


1.069 


6.43e-15 


0.80e-15 


16.73 


-3.39 


He I 


1.083 


8.94e-15 


0.80e-15 


13.2 


-4.48 


He 1+ PCygni 


1.082 


-1.95e-15 




8 


1 


Pay 


1.094 


1.27e-14 


0.16e-14 


20.5 


-6.08 


Ci(24) 


1.175 


4.70e-15 


0.90e-15 


19.8 


-2.00 


Pay? 


1.282 


2.19e-14 


O.lle-13 


21.0 


-8.53 


He I 


2.058 


1.64e-14 


0.30e-14 


33.09 


-5.163 


H2 i-os(ir 


2.122 


9e-15 


4e-15 


49.00 


-3.5 


Bry 


2.166 


2.24e-14 


0.40e-14 


25.91 


-8.87 



Line id. 


A 


F AF 


FWHM 






(fim) 


(Jyjum) 


(//m) 


(jjm) 






IRS - Spitzer 






[Sim] 


25.6 


0.04 0.008 


0.4 


-0.05 


H2 


28.2 


0.03 0.008 


0.4 


-0.02 


[Sm] 


33.5 


0.04 0.01 


0.5 


-0.05 


[Sin] 


34.8 


0.07 0.02 


0.5 


-0.07 
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Table A.13. Observed emission lines in [CTF93] 104 - ID 1 1. 



Line id. 


A 


F 


AF 


FWHM 


w,, 




()um) 


(erg cm ^ s ' ) 


A 


A 




EFOSC & Sofl - 


NTT 






Fell 


U.Ol J 


/. /Ue-1 J 


l.OOe-15 


11.3 


-0.462 


[OiJ D2- P2 


0.630 


2.00e-14 


O.lOe-14 


n 1 
9.61 


1 1 /I 
-3. 142 


[Oi] 'D2-^Pi 


0.636 


6.47e-15 


1.20e-15 


12.00 


-0.937 


Fe II 


0.646 


3.77e-15 


1.20e-15 


10.00 


-0.522 


Fe II 


0.651 


2.90e-15 


l.OOe-15 


12.00 


-0.376 


Ha 


0.656 


6.99e-13 


O.Ole-13 


11.08 


-88.21 


Hei^'Poa - 0-''Si 


0.667 


1.69e-14 


O.lle-14 


10.95 


-2.037 


[Sll] 2D5/2-%/2 + 


0.672 










[Sii] 203/2-^83/2 


0.673 


3.80e-15 


l.Oe-15 


14.1 


-0.456 


He 1 3Po2-0-3Si + 


0.706 


1.69e-14 


0.14e-14 


11.3 


-1.852 


Fe II 


0.707 










? 


0.733 


8.60e-15 


1.4e-15 


14.9 


-0.87 


? 


0.737 


4.60e-15 


1.4e-15 


14.9 


-0.468 


? 


0.743 


5.40e-15 


1.4e-15 


13.0 


-0.541 


Fe II 


0.746 


7.70e-15 


2.0e-15 


17.0 


-0.766 


Fe II 


0.771 


6.20e-14 


2.00e-14 


11.6 


-2.122 


Oi^Po-^SJ 


0.777 


3.46e-14 


0.12e-14 


11.06 


-3.504 


Pa 22 


0.836 


1 '^n^ 1 A 

1.20e-14 


0.40e-14 


14.0 


-1.324 


Pa 21 




1 A 1 A 

1.40e-14 


0.40e-14 


14.3 


-1.345 


Pa 20 


0.839 


1 Ar\^ t A 

1.49e-14 


0.40e-14 


10.1 


-1.583 


Pa 19 


0.842 


2.48e-14 


0.40e-14 


11.1 


-2.619 


Fa l6 + 


r\ QA A 
U.o44 










Oi^Po- Sj 


0.845 


/.21e-14 


0.40e-14 


1 ") 1 /I 
13.14 


-7.634 


Oi^P2-^S* 












Or^Pi-^S; 












Pa 17 


0.847 


2.63e-14 


0.40e-14 


10.8 


-2.755 


Pa 16 + 


0.850 










Cail2p3/2-2D3/2 


0.850 


2.50e-13 


0.40e-14 


10.8 


-26.5 


Cail2p3/2-2D5/2 + 


0.854 


2.64e-13 


0.40e-14 


10.8 


-27.86 


Pais 


0.855 










Pa 14 


0.860 


5.10e-14 


0.40e-14 


12.40 


-5.521 


Ni 


0.863 


9.50e-15 


0.40e-14 


15.90 


-0.984 


Can ¥1/2-^03/2 + 


0.866 


2.3e-13 


0.40e-14 


10.89 


-24.4 


Pa 13 


0.867 










Pa 12 


0.875 


6.99e-14 


0.40e-14 


11.7 


-6.23 


Mgi 


0.881 


1.29e-14 


0.40e-14 


14.6 


-1.231 


? 


0.884 


1.80e-14 


0.40e-14 


14.9 


-1.815 


Pall 


0.887 


7.00e-14 


0.40e-14 


10.9 


-7.031 


Pa 10 


0.902 


6.76e-14 


2.00e-14 


10.9 


-8.238 


Pa 9 


0.923 


8.30e-14 


2.00e-14 


12.0 


-10.2 


Pa 8 


0.955 


1.59e-13 


0.20e-13 


22.25 


-13.04 


Pa 5 


1.005 


1.94e-13 


0.20e-13 


19.45 


-14.69 


Ci(l) 


1.069 


1.70e-13 


0.08e-13 


21.20 


-10.72 


Ci(l) 


1.073 


3.44e-14 


0.60e-14 


23.2 


-1.656 


He I 


1.083 


2.88e-13 


0.07e-13 


17.6 


-17.55 


Pay 


1.094 


2.70e-13 


0.08e-13 


19.9 


-16.14 


Ci(24) 


1.175 


6.80e-14 


0.87e-14 


16.7 


-3.478 


[Pii]¥2-'D2orCi 


1.189 


3.75e-14 


0.80e-14 


18.0 


-1.90 


Si I? 


1.210 


2.59e-14 


0.80e-14 


21.8 


-1.313 


He I? 


1.270 


3.67e-14 


0.80e-14 


14.2 


-1.713 


Pay8 


1.282 


4.50e-13 


O.lle-13 


23.88 


-21.2 


Br 18 


1.535 


3.17e-14 


0.70e-14 


21.62 


-1.397 


Br 17 


1.544 


4.81e-14 


0.70e-14 


29.65 


-2.039 


Br 16 


1.556 


6.55e-14 


0.70e-14 


28.73 


-2.836 


Br 15 


1.570 


4.40e-14 


0.70e-14 


29.10 


-1.9 


Br 14 


1.588 


7.90e-14 


0.64e-14 


30.00 


-3.335 


Br 13 


1.611 


l.lle-13 


0.07e-13 


31.54 


-4.671 
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Table A.13. continued. 



Line id. 


A 


F 


AF 


FWHM 






(yum) 


(erg cm 




A 


A 


Br 12 


1.641 


1.44e-13 


0.09e-13 


36.6 


-5.878 


Ur 1 1 


1 .0(5 i 


i . / ze- 1 J 


u.uve-i J 




-O.VUj 


Nai2D3/2-5/2-'P°(m) 


1.689 


3.20e-14 


0.70e-14 


32.7 


-1.288 


Br 10 + ? 


1.737 


1.87e-13 


0.07e-13 


33.19 


-7.461 


? 


1.742 


6.23e-14 


0.70e-14 


42.5 


-2.48 


? 


L748 


4.49e-14 


0.70e-14 


25.5 


-1.88 


He I 


2.058 


1.31e-13 


0.30e-13 


33.09 


-5.163 


Bry 


2.166 


3.01e-13 


0.13e-13 


41.72 


-13.14 
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Table A.14. Observed emission lines in Meag 31 - ID 12. 



Line id. A F 


AF 


FWHM 


W,i 


(jum) (erg cm" 




A 


A 


EFOSC & Sofl 


-NTT 






H2l-0S(1) 2.122 7.2e-15 


1.5e-15 


33.0 


-4.8 


Bry 2.166 6.4e-15 


1.5e-15 


30.2 


-4.3 
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Table A.15. Observed emission lines in [CTF93] 146-2 - ID 13. 



Line id. 


A F AF 


FWHM 


w,. 




(jum) (erg cm"- s"' ) 


A 


A 




EFOSC & Sofl - NTT 






Bry" 


2.166 2.00e-15 l.OOe-15 


7 


? 



Line id. 


A 


F 


AF 


FWHM 






(jum) 


(Jyyum) 


(urn) 


(jum) 




IRS 


- Spitzer 








[Fell] 


23.3 


0.03 


0.01 


0.25 


-0.02 


H2 


28.2 


0.05 


0.01 


0.4 


-0.02 


[Sm] 


33.5 


0.15 


0.02 


0.4 


-0.04 


[Sin] 


35.3 


0.15 


0.02 


0.5 


-0.05 


[Fell] 


35.8 


0.16 


0.01 


0.4 


-0.04 


[Fei] 


36.5 


0.14 


0.07 


0.3 


-0.03 



" S/N ratio below 3cr 
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Table A.16. Observed emission lines in [CTF93] 146-1 - ID 14. 



Line id. 


A 


F 


AF 


FWHM 






(pm) 


(erg cm ^ s ') 


A 


A 




EFOSC & Sofl - 


NTT 






[Fell] aT)7/2 - (3^9/2 


1.644 


5.69e-15 


0.70e-15 


36.66 


-6.626 


[Fell] a'^Dyii-a'^Fsii" 


1.664 


1.91e-15 


0.70e-15 


33.23 


-2.074 


[Fen] a'^Dsi2 - a*Fj/2 


1.677 


3.6e-15 


0.50e-15 


39.0 


-3.8 


He I 


2.058 


1.71e-14 


0.50e-14 


46.45 


-6.946 


Bry 


2.166 


3.05e-14 


0.30e-14 


45.68 


-11.06 


[Til]? 


2.185 


6.65e-15 


2.50e-15 


32.10 


-2.375 


Nai2D3/2-5/2-'P°(m) 


2.207 


1.60e-14 


2.50e-15 


43.30 


-5.63 


CO (2-0) 


2.294 


4.06e-14 


l.OOe-14 


63.10 


-12.85 


€0(3-1)" 


2.323 


1.95e-14 


0.70e-14 


26.10 


-5.93 



Line id. 


A 


F 


AF 


FWHM 






(fim) 


(Jy//m) 


(jum) 


(jum) 




IRS 


- Spitzer 








[Sim] 


25.6 


0.09 


0.01 


0.5 


-0.03 


H2 


28.2 


0.05 


0.01 


0.4 


-0.02 


[Sm] 


33.5 


0.15 


0.02 


0.4 


-0.04 


[Sin] 


34.8 


0.1 


0.03 


0.5 


-0.3 


[Fell] 


35.8 


0.09 


0.03 


0.4 


-0.02 


[Fei] 


36.5 


0.12 


0.03 


0.4 


-0.03 



" S/N ratio below 3cr 
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Table A. 17. Observed emission lines in [CTF93] 21 1 - ID 15. 



Line id. 


A 


F 


AF 


FWHM 


w,, 




(yL(m) 


(erg cm 


-2 s-i) 


A 


A 




EFOSC & Sofl - 


IN i 1 






[Oi] D2- P2 


\J.\JD\J 


S AAp 1 s 
j.^^e-i J 


U.j(Je-i3 






[Oi] 'D2-^Pi 




z.uoe-i J 


0.50e-15 


1U.44 


-5.7 


IJ- „ 


U.U JO 


^ 1 8p 14 


r\ r\<a 1 /I 
U.Uje-i4 


11.11 


-yj.oi 


[NiJ D2- P2 




/.Vje-i J 


0.jOe-i5 


1 C 1 

l3.1 


/I 1 

-4.51 


[Sii] ^05/2-83/2 + 


U.U 11 


Ala 1 ^ 


A yt A^ 1 C 

0.40e-15 


9.8 


A yl A 

-4.49 


rc 1t\ 4c 

[Sii] D3/2- S3/2 


0.673 


3.22e-15 


A /I A^ 1 C 

0.40e-l J 


12. 1 


c c 1 
-5.51 


Do 1 1 


r\ 0/1 /I 
U.o44 










Oi^'Po-^Sj 


0.845 


2.22e-15 


0.50e-15 


9.1 


-3.4 


Ol'*P2- S* 












Oi^Pi-^S* 












Pa 16 + 


0.850 










Cail2p3/2-2D3/2 


0.850 


1.24e-14 


0.50e-15 


10.9 


-19.3 


Cail2p3/2-2D5/2-H 


0.854 


1.35e-14 


0.40e-15 


11.8 


-21.1 


Pais 


0.855 










Pa 14 


0.860 


1.86e-15 


0.60e-15 


18.80 


-2.89 


Ni 


0.863 


2.26e-15 


0.60e-15 


21.00 


-2.26 


Cail2Pi/2-'D3/2 


0.866 


1.06e-14 


0.60e-15 


10.73 


-16.5 


Pa 13 


0.867 


1.56e-15 


0.60e-15 


10.2 


-2.57 


Pa 5 


1.005 


5.52e-15 


l.OOe-15 


18.4 


-5.45 


He I 


1.083 


1.12e-14 


O.lOe-14 


12 


-8.68 


He iPCygni 


1.080 


-1.5e-15 




22.4 


1.1 


Pay 


1.094 


1.22e-14 


O.lOe-14 


30.9 


-9.46 


Pay? 


1.282 


1.97e-14 


0.20e-14 


26.14 


-10.04 


7 


1.745 


7.27e-15 


0.20e-15 


40.14 


-2.05 


H2 1-0 S(l) 


2.122 


8e-15 


2e-15 


32.29 


-2.5 


Bry 


2.166 


2.50e-14 


0.40e-15 


55.91 


-7.2 


CO? 


2.280 


9.97e-15 


3.00e-15 


20.00 


-1.9 



Line id. 


A 


F AF 


FWHM 


W,^ 




(//m) 


(Jy/im) 










IRS - Spitzer 






H2 


28.2 


0.02 0.006 


0.6 


-0.03 


[Siii]-F 


34.8 


0.1 0.03 


0.6 


-0.1 


[Siii] 


35.3 
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Table A.18. Observed emission lines in [CTF93] 187 - ID 16. 



Line id. ^ F AF FWHM 

(yum) (ergcm^^s"') A A 



EFOSC & SofI - NTT 
Bry 2.166 7.00e-14 l.OOe-14 43.2 -2.5 



Line id. 


A 


F 


AF 


FWHM 






(yum) 


(Jyjum) 


(jum) 


(jum) 






IRS ■ 


■ Spitzer 






[Sn]? 


25.2 


0.2 


0.06 


0.8 


-0.07 


[Sim] 


25.6 


0.2 


0.06 


0.6 


-0.05 


? 


26.9 


0.3 


0.06 


0.9 


-0.08 


[Sin] 


34.8 


0.13 


0.04 


0.4 


-0.04 


[Sin] 


35.3 


0.13 


0.04 


0.5 


-0.04 
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Table A.19. Observed emission lines in [CTF93] 168 - ID 17. 



Line id. 


A 


F AF 


FWHM 






(yum) 


(erg cmT^ s"') 


A 


A 






EFOSC & SofI - NTT 






Bry 


2.166 


1.80e-14 0.30e-14 


45.6 


-6.9 



Line id. 


A 


F AF 


FWHM 






(yum) 


(Jyjum) 


(jum) 


(jum) 






IRS - Spitzer 






[Sn]? 


25.2 


0.05 0.01 


0.3 


-0.05 


? 


29.7 


0.06 0.02 


0.4 


-0.01 


[Sm] 


33.5 


0.09 0.02 


0.5 


-0.08 


[Sin] 


34.8 


0.13 0.04 


0.4 


-0.04 
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Table A.20. Observed emission lines in [CTF93] 191 - ID 18. 



Line id. 


A 


F 


AF 


FWHM 


W/, 




(yum) 


(erg cm 


-2 s-i) 


A 


A 




EFOSC & Sofl - NTT 








Fe II 


0.615 


1.16e-13 


5.00e-15 


21.3 


-3.432 


Fe II 


0.617 


3.55e-14 


5.00e-15 


11.5 


-1.053 


Fe I 


0.619 


1.92e-13 


5.00e-15 


11.5 


-1.053 


Fen 


U.0Z3 


an ^ 1 A 

0.6 /e-14 


5.00e-15 


21.0 


-4.95 


[Ol] lD2-^P2 


U.o3U 


n 1 0^ 1 A 

9.1/e-14 


0.40e-14 


11.42 


-2.649 


[Oi] 'D2--^Pi 


0.636 


2.98e-14 


0.40e-14 


11.30 


-0.871 


Fe I 


0.643 


6.61e-14 


0.40e-14 


10.8 


-1.926 


Fe II 


0.646 


rvT ^1/1 

8.07e-14 


0.40e-14 


9.11 


-2.354 


Fe I 




Ci A A c 1/1 

J.44e-14 


A A C\r^ 1 A 

0.40e-14 


14.6 


-1.5a5 


Fe I 


A 1 




U.4Ue-i4 


C\ At^ 


000 
-l.wl 


[Ni] 'D2--T0 


0.653 


A A ^ 1 A 

4.43e-14 


0.50e-14 


11.8 


-1.246 


Ha 


0.656 


A m ^ 1 

4.03e-12 


0.50e-14 


9.93 


-102.9 


[Ni] 'D2--T2 


0.659 


3.84e-14 


0.50e-14 


11.4 


-0.984 


Hei ''Po2 - 0--'Si 


0.667 


3.92e-14 


0.40e-14 


10.34 


-1.043 


[Sii] 2D5/2-'*S3/2 


0.672 


1.46e-14 


0.40e-14 


14.1 


-0.395 


[Sn]2D3/2-4S3/2 


0.673 


3.04e-14 


0.40e-14 


12.3 


-0.825 


7 


0.692 


3.14e-14 


0.60e-14 


16.1 


-1.039 


Fei? 


0.697 


3.74e-14 


0.70e-14 


20.3 


-1.068 


Sii 


0.700 


2.80e-14 


0.70e-14 


13.0 


-0.845 


He I 3Po2 -0-^Si 


0.706 


3.76e-14 


0.50e-14 


13.3 


-1.107 


Fell 


0.707 


l.Ole-14 


0.50e-14 


12.2 


-0.296 


7 


0.721 


1.21e-13 


l.OOe-14 


21.9 


-4.27 


7 


0.725 


1.27e-13 


l.OOe-14 


12.4 


-4.45 


7 


0.733 


1.31e-13 


l.OOe-14 


18.4 


-4.51 


7 


0.737 


5.30e-14 


l.OOe-14 


21.6 


-1.81 


Fell 


0.746 


3.75e-14 


l.OOe-14 


11.2 


-1.248 


Fen 


0.771 


6.20e-14 


2.00e-14 


11.6 


-2.122 


Oi^Po-^SI 
Sii 


0.777 


5.98e-14 


1.50e-14 


9.53 


-2.047 


A *7A A 

0.794 


A A 1 A 

4.40e-14 


l.OOe-14 


10.7 


-1.696 


Pa 22 


O.o3o 


^ C\A ^ 1/1 

2.94e-14 


l.OOe-14 


9.8 


-1.339 


Pa 20 


A OOA 

0.839 


6.74e-14 


l.OOe-14 


10.5 


-3.03 


Pa 19 


A /f 

O.o4z 


^ AO ^ 1/1 

2.4ae-14 


l.OOe-14 


11.1 


-2.619 


Do 1 C 1 

ra io + 












Oi^Po-^S* 


0.845 


1 1 ^ IT 

1.21e-13 


l.OOe-14 


12.14 


-5.27 


Or^P2-^Sj 












Or^Pi-^S* 
Pa 17 












0.847 


6.86e-14 


l.OOe-14 


12.1 


-3.08 


Pa 16 + 


0.850 










Can2p3/2-2D3/2 


0.850 


8.70e-13 


l.OOe-14 


9.82 


-40.67 


Can2p3/2--D5/2 + 


0.854 


4.22E-13 


1.39E-15 






Pais 


0.855 










Pa 14 


0.860 


1.13e-13 


0.20e-13 


16.0 


-5.14 


Can2Pi/2-^D3/2 + 


0.866 


3.58E-13 


5.97E-15 






Pa 13 


0.867 










Pa 12 


0.875 


9.88e-14 


l.OOe-14 


9.9 


-3.07 


Mgi 


0.881 


6.32e-14 


2.00e-14 


11.5 


-3.24 


Pall 


0.887 


9.94e-14 


2.00e-14 


10.3 


-5.29 


Pa 10 


0.902 


l.OOe-13 


5.00e-14 


13.4 


-3.9 


Pa 9 


0.923 


1.60e-13 


5.00e-14 


18.0 


-10.2 


Pa 8 


0.955 


2.84e-13 


0.30e-13 


26.71 


-11.9 


Fen 


1.000 


3.25e-14 


l.OOe-14 


18.0 


-5.35 


Pa (5 


1.005 


2.92e-13 


O.lOe-13 


21.78 


-13.1 


[Sii] -P1/2-2D3/2 


1.033 


4.57e-14 


l.OOe-14 


23.1 


-2.1 


7 


1.046 


6.48e-14 


l.OOe-14 


25.4 


-2.9 


7 


1.059 


4.51e-14 


l.OOe-14 


33.1 


-2.0 


Ci(l) 


1.069 


1.81e-13 


O.lOe-13 


25.20 


-8.16 


7 


1.072 


4.58e-14 


l.OOe-14 


25.3 


-2.04 
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Table A.20. continued. 



Line id. 


A 


F 


AF 


FWHM 






(yum) 


(erg cm 


-2s-l) 


A 


A 


Ci(l) 


1.073 3.44e-14 


l.OOe-14 


23.2 


-1.656 




He I 


1.083 


2.15e-13 


A 1 (\ ,^ 1 1 

O.lOe-13 


1 T 1 

17.1 


-9.7 


He iFCygni 


10810.6 


-5.80e-14 




2.6 


1 1.6 


jti lor re ii 


1.08/ 


j.31e-14 


1 f\t\r^ 1 /I 

l.OOe-14 


OA A 

20.4 


/1 1 

-2.418 


ray 


i.uy4 


j.8e-l J 


U. lUe-1 J 


OA ^ 
/U.J 


-1 /.O 


\^ 1 ) 




7 fiflf izt 
/ .owe- 1*+ 


9 Oflf 1 /t 
z.wue- 1'+ 


91 /t 


-J. lyy 


,, 2t> 2c 

Can P3/2- 5>i/2 


1 1 0/1 

1.184 


/T/i a 1 /I 
a.o4e-14 


r\r\,3 1 /I 
/.UUe-14 


1 /.J 


A 

-4. J / 


[Pii] ■^P2- D2 or Ci 


1.189 


6.96e-14 


2.00e-14 


17.0 


-3.54 


IVfsTT 'Pi -'Do or Si T 


1.198 


6.96e-14 


2 OOe-14 


17.0 


-3.54 


Si I?" 


1.199 


7.06e-14 


2.00e-14 


23.3 


-3.6 


Si I? 


1.205 


4.82e-14 


l.OOe-14 


18.6 


-2.37 


MgiiPi-iD2 ? 


1.203 










Si I? 


1.210 


3.45e-14 


l.OOe-14 


23.3 


-1.79 


Si I? 


1.227 


3.80e-14 


l.OOe-14 


25.0 


-1.55 


Nai-P3/2--Sj^2 


1.227 










iNdi -ri/2 >Ji/2 












[i^enj a U7/2-a Ug^j 
He I + PCygni 




o OAo 1 /I 

J. sue- 14 


1 r\r\a 1 /I 
i.UUe-i4 


Zj.U 


-l.JJ 


1.270 


5.24e-14 


1 rvrV/, 1 /I 

l.OOe-14 


1 C 

15.2 


-2.88 


Pa/S 


1.282 


4.5e-13 


O.lOe-13 


24.88 


-25.2 


Mgr^P°-^^Si(m) 


1.505 


1.09e-13 


l.OOe-14 


26.6 


-7.83 


Br 20 


1.519 


2.30e-14 


l.OOe-14 


30.9 


-1.67 


Br 17 


1.544 


2.06e-14 


l.OOe-14 


34.24 


-1.55 


Br 16 


1.556 


1.64e-14 


l.OOe-14 


22.24 


-1.25 


Nai^Si/2-^P"(m) Mgi 


1.576 


7.53e-14 


l.OOe-14 


60.6 


-5.72 


Br 14 


1.588 


5.08e-14 


l.OOe-14 


35.50 


-3.84 


9 


1.600 


3.40e-14 


l.OOe-14 


32.10 


-2.6 


r>r i J 


1 1 1 
l.ol i 


0. /ue-14 


1 C\C\a 1 A 

i.uue-i4 


T/i ^/i 


A f^n 1 
-4.0 / 1 


Br 12 


1.641 


4.21e-14 


O.lOe-13 


27.9 


-3.2 


[Fell] a'*D7/2-a'*F9/2 


1.644 


3.81e-14 


O.lOe-13 


33.7 


-3.2 


7 


1.667 


2.53e-14 


0.40e-14 


31.20 


-1.966 


Br 1 1 


1.681 


6.19e-14 


0.40e-14 


35.61 


-5.93 


_/ 


1 /CO yl 


/I T 1/1 
/.43e-14 


A A A.^ 1 A 

0.40e-14 


"3 1 AA 
31.00 


1 A 1 

-1.91 


i>ai U3/2-5/2 r I^IIi; .' 


i .VJO" 










BrlO+ ? 


1.737 


1.26e-13 


0.07e-13 


52.29 


-10.65 


7 


1.742 


2.78e-14 


0.70e-14 


28.5 


-2.38 


He I 


2.058 


9.07e-14 


0.30e-13 


35.23 


-10.38 


Br 7 


2.166 


1.41e-13 


0.05e-13 


36.8 


-16.54 


Nai2D3/2-5/2-'P°(m) 


2.202 


2.53e-14 


0.50e-14 


42.16 


-3.29 


Nai2D3/2-5/2-'P"(m) 


2.207 


4.07e-14 


0.50e-14 


50.00 


-5.309 



Line id. 


A 


F 


AF 


FWHM 






(yum) 


(Jy//m) 




{^im) 




IRS- 


Spitzer 








[Sim] 


25.6 


0.06 


0.01 


0.5 


-0.03 




28.2 


0.028 


0.01 


0.3 


-0.01 


[Sm] 


33.5 


0.09 


0.02 


0.4 


-0.04 


[Sin] 


34.8 


0.07 


0.02 


0.4 


-0.03 


[Fen] 


35.8 


0.1 


0.02 


0.4 


-0.04 
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Table A.21. Observed emission lines in [CTF93] 246 A - ID 19. 



Line id. 


A 


F AF 


FWHM 






(//m) 


(erg cm"- s"') 


A 


A 






EFOSC & SofI - NTT 






Bry 


2.166 


3.0e-15 l.OOe-15 


46.3 


-9.8 



Line id. 


A 


F AF 


FWHM 


w,, 




(yum) 


(Jyjum) 


(jum) 


(fim) 






IRS - Spitzer 






[Sim] 


25.6 


0.05 0.01 


0.4 


-0.12 


[Sm] 


33.5 


0.04 0.01 


0.5 


-0.08 


[Fei] 


36.5 


0.06 0.02 


0.4 


-0.13 
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Table A.22. Observed emission lines in [CTF93] 186 - ID 20. 



Line id. 


A 
(yum) 


F AF 
(erg cmT^ s"') 


FWHM 
A 


A 


Bry 


2.166 


EFOSC & SofI - NTT 
1.90e-14 0.30e-14 


54.2 


-7.1 




Line id. 


A 
(yum) 


F AF 
(Jyjum) 


FWHM 

(jum) 


w, 

(jum) 


[Siii] 


33.5 


IRS - Spitzer 
0.15 0.02 


0.5 


-0.05 
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Table A.23. Observed emission lines in [CTF93] 246 B - ID 21. 



Line id. 


A 
(//m) 


F AF 
(erg cm"- s"') 


FWHM 
A 


A 


Bry 


2.166 


EFOSC & SofI - NTT 
7.5e-15 l.OOe-15 


51.7 


-7.3 




Line id. 


A 
(yum) 


F AF 
(Jyjum) 


FWHM 

(jum) 


w,, 

(urn) 


[Sim] 
[Sin] 


25.6 
34.8 


IRS - Spitzer 
0.015 0.02 
0.09 0.02 


0.3 
0.3 


-0.03 
-0.2 
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Table A.24. Observed emission lines in [CTF93] 237-2 - ID 22. 



Line id. 


A F AF 
(/urn) (erg cm"- s"') 


FWHM 
A 


A 




EFOSC & SofI - NTT 






H2 1-0 S(l) 


2.122 2.0e-15 l.Oe-15 


38.0 


-2.8 


Bry" 


2.166 <1.0e-15 







" S/N ratio below 3cr 
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Table A.25. Observed emission lines in [CTF93] 237-1 - ID 23. 



Line id. /i F AF 


FWHM 




(yum) (erg cm"^ s"') 


A 


A 


EFOSC & Sofl - NTT 






H2l-0S(1) 2.122 8.1e-15 2.5e-15 


38.0 


-3.4 


Bry 2.166 1.20e-14 0.30e-14 


29.0 


-5.0 
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Table A.26. Observed emission lines in [CTF93] 216-1 - ID 24. 



Line id. 


A 


F 


AF 


FWHM 






(yum) 


(erg cm 




A 


A 




EFOSC & Sofl - NTT 






L^enj a U^ii - at^ii 


1 


7 OQp-l S 


W.UWC- i »7 


J J. i i 


-i J.J 


yte iij a Z^3/2 - a t-ii2 


1 fiflO 
i .DUW 


o.uue-i J 


W.oUc- i J 


11 1 1 
J J. i i 


-i J.J 


[Fe ii] a^D-112 — 0^/^9/2 




9 zt8f> 1 zt 


A Kc\(:* 1 

U. JUc-i J 




zt«; s 

-^j.o 


[re rrj a /J1/2 — a r 


1.0o4 


J. j4e-l J 


U.4Ue-i J 


9T 1 9 

Z/ AZ 


-o.Oj 


[Fe 11] a^D5/2 - 0^^7/2 


1.0/7 


*7 /I 1 C 

7.46e-15 


0.50e-15 


33.01 


-13.2 


[Fe 11] flX)3/2 - 13 ^5/2 


1 *7 1 

1.712 


A 1 ^ 1 C 

3.01e-15 


A CA^ 1 C 

0.50e-15 


35.21 


-4.89 


H2 1-0 !S(o) 


1 -7 1/1 

1. /14 


1 f\r\^ 1 c 

l.OOe-15 


0.40e-l J 


1 T AA 

17.00 


A OA 

-4.89 


[Fe 11] a^Pi - a^F4 ? 


1. /jy 


1 If^a 1 ^ 

l.joe-lj 


r\ A r\a 1 ^ 

U.4Ue-l J 


31.60 


-1.81 


[Feu] a'^Di/2 - a'^F3/2 + 


1.745 










H2 1-0 S(7) 


1.749 


1.51e-15 


0.40e-15 


27.30 


-1.95 


[Fell] a'^Ds/2- a'^Fspf 


1.800 


5.36e-15 


2.00e-15 


30.60 


-5.88 


TT 1 r\C/'1\ 

H2 1-0 S(l) 


2.122 


6.86e-15 


0.60e-15 


T A A A 

29.40 


A AA 

-4.90 


[Fell] a^P3/2-a%/2 


2.133 


2.11e-15 


0.60e-15 


33.90 


-1.5 


Biy 


2.166 


1.74e-14 


0.80e-15 


49.06 


-12.42 


JTii] 
[Fen] a^Hcii2 - a^7/2 


2.179 


3.20e-15 


0.80e-15 


29.2 


-2.26 


2.254 


7.02e-I5 


1.50e-15 


38.59 


-4.6 


CO (2-0) 


2.294 


2.15e-14 


0.50e-14 


52.11 


-12.92 


CO (3-1) 


2.323 


2.13e-14 


0.50e-14 


37.83 


-12.2 




IRS - Spitzer 








[Sm] 


33.5 


0.16 


0.02 


0.4 


-0.02 


[Sin] 


35.3 


0.13 


0.03 


0.3 


-0.2 



" S/N ratio below 3cr 
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Table A.27. Observed emission lines in [CTF93] 216-2 - ID 25. 



Line id. 


A F AF 


FWHM 






(jum) (erg cm"^ s"') 


A 


A 




EFOSC & Sofl - NTT 






Bry 


2.166 9.6e-14 l.Oe-14 


59.05 


-3.6 
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Table A.28. Observed emission lines in [CTF93] 245 B-2 - ID 26. 



Line id. 


A 


F 


AF 


FWHM 






(yum) 


(erg cm 




A 


A 




EFOSC & Sofl - NTT 






[Fen] ai)7/2 - a /^9/2 


1.644 


2.32e-14 


0.50e-15 


T r\f\ 

28.00 


At T 1 

-41.21 


[Fell] flT)i/2 - 0^5/2 


1.664 


3.58e-14 


0.50e-15 


28.93 


-6.50 


[Fe n] a^Ds/2 - a^Fy2 


1.677 


6.67e-15 


0.50e-15 


29.87 


-11.59 


[re iij a IJ^/') ~ a t 


1.712 


2.94e-15 


0.50e-15 


3o. 1 / 


-4.00J 


H2 1-0 8(8)" 


1.714 


l.OOe-15 


0.40e-15 


17.00 


-4.685 


[Fell] - a'^F^/2 + 


1.745 










H2 1-0 S(7) 


1.749 


1.68e-15 


0.50e-15 


30.90 


-2.166 


H2 1-0 S(l) 


2.122 


6.35e-15 


0.60e-15 


26.60 


-4.51 


[Fell] fl2p3/2 - fl'*P3/2 


2.133 


2.12e-15 


0.60e-15 


36.72 


-1.498 




2.166 


1.56e-14 


0.80e-14 


44.72 


-11.06 



" S/N ratio below 3cr 
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Table A.29. Observed emission lines in [CTF93] 245 B-1 - ID 27. 



Line id. 


A 


F 


AF 


FWHM 






(yum) 


(erg cm 




A 


A 






EFOSC & SofI - NTT 






Pad 


1.005 


2.31e-15 


0.80e-15 


24.25 


-3.78 


He I 


1.083 


3.96e-15 


0.90e-15 


20.5 


-3.94 


Pay 


1.094 


3.13e-15 


l.OOe-15 


24.1 


-2.96 


? 


1.135 


4.56e-15 


l.OOe-15 


23.17 


-3.37 


Ci? 


L166 


4.25e-15 


l.OOe-15 


28.19 


-2.76 


Ci 


L175 


2.97e-15 


l.OOe-15 


20.0 


-1.88 


? 


1.278 


8.24e-15 


2.00e-15 


29.5 


-3.3 


Pa/3 


1.282 


1.43e-14 


2.00e-15 


28.33 


-5.5 


Bry 


2.166 


4.37e-14 


0.60e-14 


58.97 


-5.2 
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o 
o 
o 



u 
Q 




RA (J2000) 



Fig. B.l. Spitzer/IRAC three-colour image (3.6, 4.5, and S.Oyum - blue, 
red, and green, respectively) of the [CTF93]216 system, sources #24 
(-1), and #25 (-2), along with their jets (in green). These images were 
taken in 2006, before the outburst of source #25. 



Appendix B: Spitzer-IRAC flows in [CTFg3]216 

Figure B.l shows the Spitzer/IRAC three-colour image (3.6, 
4.5, and S.Ojum - blue, red, and green, respectively) of the 
[CTF93]216 binary flow system. The two precessing jets from 
source #24 ([CTF93] 216-1), and #25 ([CTF93] 216-2) are vis- 
ible in green, and they are roughly E-W, and NE-SW oriented, 
respectively. 



